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Introduction 

As  indicated  in  the  previous  annual  progress  report,  limited  success  was  found  employing 
the  fuzzy  clustering  algorithm  on  the  first  43  subjects  studied.  Clustering  temporal  pixel 
image  intensity  curves  based  on  a  series  of  four  to  six  time  points  post  gadolinium 
injection  presented  a  limited  set  of  functional  variants.  A  more  direct  physiological 
approach  toward  understanding  sodium  content  in  breast  lesions,  which  is  one 
physiological  parameter  that  regulates  temporal  gadolinium  uptake.  Special  coils  and 
pulse  sequences  were  developed  to  perform  sodium  breast  imaging  on  a  clinical  MRI 
scanner.  The  resulting  sodium  breast  images,  included  in  this  document  were  presented  at 
the  recent  Era  of  Hope  meeting.  The  factors  that  influence  sodium  content  in  a  benign  or 
malignant  lesion  derive  from  several  sources.  Palpable  lesions  likely  derive  stiffness 
(relative  to  surrounding  tissue)  through  osmotic  pressure  buildup  in  the  local  extra¬ 
cellular  matrix.  The  most  studied  tissue  where  sodium  concentration  is  measured  in  this 
regard  is  articular  cartilage.  Thus,  cartilage  has  been  employed  as  a  convenient  model 
tissue  for  quantifying  the  osmotic  mechanisms  associated  with  elevated  sodium  content. 

Background 

Studies  have  indicated  an  interstitial  fluid  pressure  of  29  mm  Hg  in  invasive  ductal 
carcinomas,  3.6  mm  Hg  in  benign  lesions,  and  -0.3  mm  Hg  in  normal  breast  tissue  (P.  B. 
Roemer  et.  al.,  MRM  16:192,1990).  An  increase  in  biological  fluid  pressure  is  often 
osmotically  related  to  an  increase  in  tissue  sodium  content.  To  non-invasively  explore  the 
notion  that  high  sodium  content  is  present  in  malignant  breast  lesions,  technology  for 
improving  the  sensitivity  and  resolution  of  sodium  breast  MR  images  must  be  developed. 
We  have  developed  a  transmit  coil/  phased  array  system  for  sodium  to  be  used  at  4  Tesla 
By  double  tuning  the  transmit  coil  and  building  special  phased  arrays  with  both  proton 
and  sodium  coils,  we  perform  breast  imaging  at  both  proton  and  sodium  frequencies 
within  the  same  examination. 

Magnetic  resonance  imaging  of  breast  lesions  has  been  unable  to  identify  malignant  vs. 
benign  lesions  using  only  T1 -weighted  and  T2- weighted  imaging  sequences.  Presently, 
gadolinium-based  contrast  agents  are  being  explored  as  an  approach  toward 
distinguishing  lesion  types.  The  increased  accumulation  of  contrast  agent  in  malignant 
tumors  is  thought  to  arise  from  the  increased  permeability  of  the  blood  vessels  within  the 
tumor.  This  increased  permeability  combined  with  the  lack  of  a  lymphatic  system  may 
also  contribute  to  an  increase  in  interstitial  fluid  pressure  and  an  increase  in  tissue  sodium 
content. 


Coil  Design 

The  transmit  coil  was  made  with  0.6  cm  wide  copper  tape  in  a  biplanar  design.  The 
geometry  is  shown  in  the  figure.  This  coil  makes  linear  Bl.  The  breast  falls  into  the  coil 
from  die  top  and  the  receiver  coils  are  used  to  lightly  compress  the  breast  from  the  sides. 
The  capacitance  for  the  transmit  coil  is  distributed  to  six  sites.  The  coil  is  driven 
capacitively.  In  order  to  detune  the  transmit  coil  during  receive  conditions,  a  diode  is 
placed  in  series  with  the  drive  of  the  coil.  This  diode  is  biased  on  during  transmit 


conditions  with  a  5V  dc  signal  which  is  delivered  by  the  spectrometer  with  the  rf  pulse. 
It  is  back  biased  with  -  15V  dc  during  all  other  times. 


Transmit  Coil  Geometry 


Two  plates  each  with  two  receive  coils  were  used  to  lightly  compress  die  breast  during 
the  examination  in  die  four  coil  array  set  up.  The  dimensions  are  shown  in  the  receiver 
coil  geometry  figure.  The  coils  are  etched  on  single  side  copper  clad  epoxy  circuit  board. 
The  capacitance  was  distributed  at  6  locations.  A  detuning  circuit  as  described  by 
Roemer,  et  al,  was  used  for  decoupling  the  receiver  coils  from  die  transmitter  coil  during 
the  excitation  pulse. 


Receiver  Coil  Geometry 


Results  and  Discussion 

Figures  1  and  2  show  an  axial  proton  and  sodium  image  of  a  normal  breast  Even  though 
there  are  no  abnormalities  in  this  breast,  variation  in  the  sodium  signal  can  be  seen.  For 
instance,  the  regions  of  lipid  signal  do  not  show  much  sodium  signal  while  the  breast 
parenchyma  does.  In  die  sagittal  image  (figure  3),  a  gap  in  the  sodium  signal 
representative  of  the  bones  in  the  chest  can  be  seen.  In  addition,  variations  in  the  level  of 
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the  sodium  signal  is  evident  which  is  indicative  of  the  variations  in  the  parenchymal 
thickness. 

In  normal  breast  tissue,  sodium  signal  can  be  seen  in  parenchymal  tissue  with  variations 
most  likely  indicative  of  parenchymal  density.  Skin  also  produces  a  signal  as  seen  on  the 
axial  images.  The  sodium  images  need  to  be  acquired  at  fairly  high  spatial  resolution  (1- 
7mm)  to  maintain  correspondence  with  the  proton  images,  which  have  an  in-plane 
resolution  of  0.4mm.  This  correspondence  will  permit  the  comparison  of  sodium  content 
with  contrast  agent  uptake. 
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Figure  1.  Axial  proton  image  of  human  breast.  The 
proton  image  is  a  rf- spoiled  gradient  echo  with  a  FOV  of 
24  cm  (only  part  of  the  FOV  shown  here),  TE/TR  = 
2.5/100  and  7  mm  slice  thickness. 


Figure  2.  Axial  sodium  image  of  human  breast.  The 
sodium  image  was  acquired  fast,  rf  spoiled  gradient 
echo  sequence  was  used  with  a  flip  of  90  degrees, 
TE/TR  =  2.5/150  ms,  15.6  kHz  receiver  bandwidth, 
acquisition  matrix  of  256x128,  sixteen  summations  and 
a  FOV  of  30  cm  and  slice  thickness  of  30  mm. 


Figure  3.  Sagittal  sodium  image  of  human  breast.  This  was  acquired  with  a 
bandwidth  of  7.81  kHz  and  8  summations  resulting  in  an  acquisition  time  of  2  min. 
35  sec.  A  fast,  r£  spoiled  gradient  echo  sequence  was  used  with  a  flip  of  90 
degrees,  TE/TR  =  3.3/150  ms,  acquisition  matrix  of  256x128,  a  FOV  of  30  cm  and 
slice  thickness  of  30  mm. 


KEY  RESEARCH  ACCOMPLISHMENTS 


•  Realized  sodium  MRI  of  the  human  breast. 

•  Modeled  the  contribution  of  sodium  equilibria  to  the  physical 
properties  of  connective  tissue. 
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REPORTABLE  OUTCOMES 


The  study  of  sodium  imaging  of  the  human  breast  was  presented  at  the  Era  of  Hope 
Meeting.  A  study  of  sodium  equilibria  in  connective  tissue  was  published  as  a 
dissertation. 

•  Sodium  MRI  of  the  human  breast,  Kaufman  JH,  Bolinger  L,  Leigh  JS, 
University  of  Pennsylvania  Medical  Center,  presented  at  US  Army  Era  of 
Hope  Meeting,  Atlanta  GA,  2000. 

•  Kaufman  JH,  Ph.D.  Dissertation,  Biochemistry  and  Molecular  Biophysics, 
University  of  Pennsylvania  School  of  Medicine.  December  2000. 
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SODIUM  MRI  OF  THE  HUMAN  BREAST 


Kaufman  JH,  Bolinger  L,  Leigh  JS 


University  of  Pennsylvania  Medical  Center 


jonathak@mail.med.upenn.edu 


Breast  tumors  have  a  higher  sodium  content  than  normal  tissue.  Thus,  sodium  MRI  could 
prove  to  be  a  valuable  non-invasive  diagnostic  for  discriminating  malignant  lesions.  We 
have  developed  a  technique  for  generating  sodium  images  of  the  breast  in  a  clinical 
setting. 

A  sodium  transmitter  with  mutlicoil  receive  was  built  of  a  similar  design  to  the  one  used 
for  proton  breast  examinations  at  4  Tesla.  The  transmitter  is  a  biplanar  design,  which  is 
easily  double-tuned  to  both  proton  and  sodium.  Two  plates  each  with  two  receive  coils 
were  used  to  lightly  compress  the  breast  during  the  examination. 

Figure  1  shows  an  axial  proton  and  sodium  image  of  a  normal  breast.  Even  though  there 
are  no  abnormalities  in  this  breast,  variation  in  the  sodium  signal  can  be  seen.  For 
instance,  the  regions  of  lipid  signal  do  not  show  much  sodium  signal  while  the  breast 
parenchyma  does.  In  the  sagittal  image  (figure  2),  a  gap  in  the  sodium  signal 
representative  of  the  bones  in  the  chest  can  be  seen.  In  addition,  variations  in  the  level  of 
the  sodium  signal  is  evident,  which  is  indicative  of  the  variations  in  the  parenchymal 
thickness. 


Figure  1(a)  Figure  1(b)  Figure  2 


The  U.S.  Army  Medical  Research  and  Material  Command  under  DAMD17-98-1-8231 
supported  this  work. 
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Abstract 


Magnetic  Resonance  Studies  of  Cartilage  Osmotic  and  Mechanical 

Properties 

Jonathan  H.  Kaufman 
John  S.  Leigh 

Articular  cartilage,  the  tissue  that  coats  the  ends  of  bones  in  synovial  joints,  is  remarkable 
in  that  it  is  solid  while  being  primarily  composed  of  water.  The  progressive  deterioration 
of  articular  cartilage  debilitates  over  20%  of  the  human  population.  Most  major  pharma¬ 
ceutical  companies  actively  pursue  cartilage  regenerative  therapies;  however,  it  is  diffi¬ 
cult  to  test  the  effectiveness  of  novel  treatments  since  cartilage  condition  cannot  be  moni¬ 
tored  non-invasively.  Here,  we  present  several  new  ways  to  assess  cartilage  softening  via 
magnetic  resonance  imaging  (MRI).  Cartilage  softening  is  one  of  the  earliest  manifesta¬ 
tions  of  cartilage  degeneration.  First,  a  new  quantitative  osmotic  model  of  cartilage 
rheology  is  developed.  Then,  the  equilibrium  of  a  charged  MR  contrast  agent  is  used  to 
confirm  the  osmotic  properties  that  are  hypothesized  to  give  cartilage  its  strong  compres¬ 
sive  resistance.  A  novel  method  for  mechanically  compressing  cartilage  specimens  in  an 
imaging  environment  is  then  demonstrated.  Conventional  MRI  parameters  of  articular 
cartilage  are  presented  as  a  function  of  compression.  One  major  advantage  that  this  new 
technique  has  over  purely  mechanical  studies  is  that,  during  compression,  positional 
cartilage  water  content  can  be  dynamically  mapped.  Additional  techniques  are  presented: 
one  for  dynamically  measuring  cartilage  fluid  permeability,  and  a  second  for  spatially 
mapping  cartilage  rheological  properties.  The  osmotic  model  is  then  further  confirmed  by 
interleaved  sodium  and  proton  imaging  dining  compression. 


v 


Contents 


Chapter  1  A  review  of  articular  cartilage  structure,  function,  and 
imaging 

1 

Chapter  2  A  magnetic  resonance  imaging  based  demonstration  of 
articular  cartilage  ionic  equilibria 

10 

Chapter  3  A  method  for  compressing  articular  cartilage  within  a 
magnetic  resonance  imaging  environment 

29 

Chapter  4  Spectroscopic  magnetic  resonance  imaging  parameters  of 
articular  cartilage  during  compression 

38 

Chapter  5  The  cartilage  uni-axial  confined  compression  experiment 
within  a  magnetic  resonance  imaging  environment 

62 

Chapter  6  Articular  cartilage  fluid  permeability  during  compression 

72 

Chapter  7  Multinuclear  study  of  the  cartilage  osmotic  model 

79 

Appendix  A  Magnetic  resonance  imaging  of  unconfined  cartilage 
compression 

97 

Appendix  B  Magnetic  resonance  imaging  cartilage  strain  mapping 
during  compression 

108 

Bibliography 

139" 

v 


List  of  Tables 


Table  2.1  Ti  measurements  of  the  five  cartilage  specimens  equilibrated 
in  137  mM  PBS  without  Gd(DTPA). 

22 

Table  2.2  Ti  measurements  of  the  five  cartilage  specimens  equilibrated 
in  PBS  solutions  of  various  NaCl  concentrations  with  1.0  mM 

Gd(DTPA). 

25 

Table  4.1  Proton  relaxation  parameters.  Relaxation  times  are  given  for 
proton  Ti,  and  proton  T2  for  cartilage  in  the  normal  and  mechanically 
compressed  state. 

51 

Table  4.2  The  statistical  significance  of  the  change  in  relaxation  time 
constants  associated  with  compression. 

52 

Table  4.3  Proton  MR  mean  relaxation  time  constants  of  normal  and 
proteoglycan  depleted  cartilage  specimens. 

61 

Table  5.1  The  deformation  creep  data  of  figure  6.3  were  fit  to  the 
biexponential  decay  function  given  in  equation  6.2. 

70 

Table  7.1  Sodium  and  proton  relaxation  parameters. 

91 

Table  7.2  The  statistical  significance  of  the  change  in  relaxation  time 
constants  associated  with  compression. 

92 

vi 


Chapter  1 


Articular  Cartilage  Structure,  Mechanics  and  Imaging 

1.1  Overview 

This  thesis  describes  several  magnetic  resonance  imaging  (MRI)  experiments  per¬ 
formed  on  articular  cartilage.  Each  experiment  was  designed  to  elucidate  articular  carti¬ 
lage  rheological  mechanisms,  especially  with  regard  to  changes  associated  with  os¬ 
teoarthritis.  The  intention  of  this  work  is  to  propose  and  test  a  new  theoretical  model  for 
cartilage  deformation.  The  primary  need  for  a  new  model  stems  from  the  inadequacy  of 
existing  cartilage  deformation  models,  which  are  either  inaccurate  or  overly  complex. 

This  section  begins  with  a  brief  description  of  articular  cartilage  anatomy  and  bio¬ 
chemistry  in  health  and  disease,  specifically  with  regard  to  osteoarthritis.  Next,  a  working 
hypothesis  for  articular  cartilage  deformation  is  presented  along  with  a  strategy  to  explore 
the  theoretical  implication  of  the  hypothesis  as  well  as  steps  toward  experimentally  test¬ 
ing  the  consequential  model.  The  new  model  is  presented  in  light  of  the  most  commonly 
accepted  existing  model  for  cartilage  deformation.  An  overview  of  the  thesis  progression 
is  then  provided. 

1.2  Articular  Cartilage 

The  primary  function  of  articular  cartilage  is  to  provide  a  low  friction  surface  be¬ 
tween  bones  in  synovial  joints.  Articular  cartilage  is  composed  of  water,  a  network  of 
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collagen  fibrils,  and  proteoglycan  aggregates.  Compared  to  most  tissues,  cartilage  is  rela¬ 
tively  acellular  (McCutchin  et  al.  1972). 

Approximately,  80%  of  cartilage  (by  weight)  is  water;  the  remainder  constitutes 
the  solid  matrix.  The  solid  matrix  is  primarily  composed  of  two  biopolymers  (type-II  col¬ 
lagen  fibrils,  and  proteoglycan  aggregates)  which  serve  distinct  functions.  Having  a  high 
tensile  strength,  the  collagen  fibrils  maintain  cartilage  integrity.  The  proteoglycan  aggre¬ 
gates  impart  a  substantial  negative  fixed  charge  density  within  the  cartilage  matrix  which 
causes  a  large  osmotic  pressure.  The  counteracting  forces  of  the  high  osmotic  pressure 
and  the  collagen  fibril  tensile  strength  give  cartilage  its  solid  character  (Mankin  et  al. 
1994).  This  concept  is  discussed  in  more  detail  and  quantified  in  chapter  2. 

The  high  water  pressure  within  cartilage  presents  a  resistance  to  mechanical  com¬ 
pression  and  permits  articular  cartilage  to  serve  as  a  mechanical  shock  absorber  when 
synovial  joints  undergo  brief  periods  of  impact.  During  cartilage  compression,  cartilage 
fluid  is  released  which  contributes  to  synovial  joint  friction  reduction. 

Articular  cartilage  is  the  form  of  cartilage  that  resides  at  the  articulating  surfaces 
of  bones  in  synovial  joints.  Several  other  types  of  cartilage  exist.  For  example  one  form 
of  cartilage  that  also  exists  in  certain  synovial  joints  is  the  cartilage  that  form  the  meiscus 
in  the  human  knee  joint.  Another  form  of  cartilage  involved  in  joint  mobility  is  the  carti¬ 
lage  that  forms  the  discs  in  a  vertebral  column.  A  more  structurally  distinct  type  of  carti¬ 
lage  that  exists  outside  of  synovial  joints  is  structural  cartilage.  This  type  of  cartilage  is  a 
major  component  of  certain  solid  structures  in  the  human  body  such  as  the  trachea,  nose 
and  external  ear. 
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One  unique  aspect  of  articular  cartilage  is  its  layered  anatomy.  Articular  cartilage 
is  often  described  in  the  histological  literature  as  having  three  relatively  distinct  layers. 
The  layer  at  the  articulating  surface  is  usually  designated  the  superficial  zone.  In  this 
zone,  the  collagen  fibers  typically  run  parallel  to  the  articular  surface.  The  layer  adjacent 
to  the  subchondral  bone  is  usually  designated  the  radial  zone.  In  this  zone,  collagen  fibers 
are  typically  oriented  perpendicular  to  the  cartilage-bone  interface.  The  radial  zone  is 
usually  highly  calcified.  The  region  of  articular  cartilage  between  the  superficial  zone  and 
the  radial  zone  is  usually  referred  to  as  the  intermediate  zone.  In  this  zone,  collagen  fibers 
do  not  have  a  distinct  orientation.  Scanning  electron  microscopy  studies  of  articular  carti¬ 
lage  yield  an  overall  sheet-like  structure.  In  these  studies,  cartilage  structure  is  in  the  form 
of  parallel  sheets  that  are  oriented  parallel  to  the  articular  surface  in  the  superficial  zone, 
and  aggregately  bend  in  a  direction  perpendicular  to  the  surface.  (Mow  et.  al.  1990). 

The  notion  of  a  layered  cartilage  structure,  and  the  notion  that  cartilage  structure 
is  locally  non-isotropic,  is  important  to  experiments  and  discussions  described  in  the  pre¬ 
sent  work.  Specifically,  the  present  work  employs  the  assumption  that  cartilage  has  local 
lateral  (parallel  to  the  articular  surface)  mechanical  homogeneity,  but  has  mechanical 
non-homogeneity  in  a  direction  perpendicular  to  the  articular  surface. 

1.2  Osteoarthritis 

Osteoarthritis  is  a  condition  present  in  most  (if  not  all)  species  that  contain  syno¬ 
vial  joints  (Lawrence  et  al.  1989).  The  disease  is  characterized  by  a  progressive  deteriora¬ 
tion  of  articular  cartilage.  There  is  evidence  that  the  early  onset  of  osteoarthritis  is  charac¬ 
terized  by  the  selective  loss  of  the  proteoglycan  component  of  the  solid  matrix.  There  is 
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currently  no  pharmacological  treatment  for  reversing  the  progression  of  osteoarthritis 
(Ghosh  et  al.  1991).  Existing  therapies  include:  pain  management  (Davis  et  al.  1993, 
Todd  et  al.  1990),  joint  strengthening  exercises  (Jan  et  al.1993,  Kovar  et  al.  1992),  dietary 
supplements,  prosthetic  joint  replacement  (Amendolaet  al.  1989),  and  chondrocyte  trans¬ 
plantation.  One  major  impediment  to  pharmaceutical  development  for  therapies  that  re¬ 
verse  the  osteoarthritic  process  is  that  the  progression  or  regression  of  osteoarthritis  is 
currently  not  measurable  non-invasively  (March  et  al.1996).  It  is  therefore  difficult  to 
measure  the  efficacy  of  novel  therapeutics. 

The  disease  process  of  osteoarthritis  involves  the  degeneration  of  articular  carti¬ 
lage.  The  earliest  putative  changes  in  cartilage  due  to  osteoarthritis  are  softening  of  carti¬ 
lage,  a  decrease  in  proteoglycan  content,  and  an  increase  in  permeability  (Armstrong  et 
al.  1982;  Lesperance  et  al.  1992;  Mansour  et  al.  1976).  Each  of  these  changes  has  been 
associated  with  an  alteration  in  the  mechanical  properties  of  cartilage  (Mow  et  al.  1984; 
Basser  et  al.  1988;  Mow  et  al.  1980).  Therefore,  there  has  been  substantial  interest  in 
studying  the  mechanical  properties  of  normal  and  degenerated  articular  cartilage. 

The  impact  and  motivation  of  the  efforts  described  in  this  thesis  stem  from  the 
substantial  prevalence  of  osteoarthritis.  Approximately  40  million  people  (in  the  U.S. 
alone)  are  affected  by  osteoarthritis,  with  a  major  characteristic  being  a  change  in  the  me¬ 
chanical  properties  of  articular  cartilage.  Since  there  is  currently  no  method  for  diagnos¬ 
ing  osteoarthritis  in  its  early  stage,  as  well  as  no  current  method  for  sensitively  measuring 
the  ongoing  progression  of  osteoarthritis  in  any  given  individual,  a  new  diagnostic 
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method  is  needed.  An  effective  approach  may  likely  include  consideration  of  cartilage 
mechanical  properties. 

1.3  Cartilage  Mechanical  Properties 

Early  attempts  to  study  the  mechanical  properties  of  articular  cartilage  involved 
indentation  experiments.  In  these  experiments,  the  dynamics  of  cartilage  response  to  an 
indenter  was  used  to  evaluate  cartilage  compressive  modulus  and  cartilage  viscosity 
(Harris  et  al  1972).  The  quantitative  results  of  such  experiments  were  difficult  to  interpret 
since  cartilage  is  non-isotropic  and  indentation  experiments  are  inherently  multidimen¬ 
sional.  The  results  of  these  experiments  varied  when  different  shaped  indenters  were 
used.  Also,  the  shape  of  an  indentation  at  any  given  time  was  difficult  to  characterize. 
One  would  intuitively  expect  that  observing  the  response  of  cartilage  to  indentation 
would  yield  information  on  its  mechanical  properties.  However,  in  practice,  indentation 
experiments  have  been  unable  to  yield  substantial  quantitative  information  on  cartilage 
mechanics.  The  primary  reason  for  this  lack  of  effectiveness  is  the  complexity  involved 
with  developing  a  theory  to  adequately  describe  cartilage  response  to  indentation,  espe¬ 
cially  given  the  fact  that  cartilage  response  to  indentation  depends  on  the  shape  of  the  in- 
dentor. 

This  challenge  was  met  by  the  design  of  the  uniaxial  confined  compression  ex¬ 
periment  (Mow  et  al.  1980).  Even  though  indentation  experiments  have  maintained  utility 
for  testing  certain  cartilage  deformation  hypotheses  (Mow  et  al.  1989;  Jurvelin  et  al. 
1990;  Oloyede  et  al.  1991;  Hale  et  al.  1993;  Athanasiou  et  al.  1994;  Obeid  et  al.  1994; 
and  Haut  et  al.  1995),  the  advent  of  the  uniaxial  compression  experiment  yielded  experi- 
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mental  results  that  could  be  directly  related  to  a  theoretical  description  of  cartilage  uniax¬ 
ial  compression.  It  should  be  noted  that  the  uniaxial  compression  experiment,  of  course, 
did  not  solve  all  challenges.  For  example,  it  has  also  been  shown  that  the  dynamics  of 
cartilage  compression  are  dependant  upon  load  cycling  (Suh  et  al.  1995;  Su  et  al.  1996). 
Moreover,  the  theoretical  cartilage  mechanical  model  for  the  uniaxial  compression  ex¬ 
periment  (as  well  as  other  experiments)  was  still  fairly  complex. 

In  addition  to  compression  studies,  several  other  types  of  mechanical  studies  have 
been  used  to  evaluate  cartilage  mechanical  properties.  Three  common  methods  are  steady 
state  permeation,  stress  relaxation,  and  torsional  shear  (Zahler  1999).  In  the  steady  state 
permeation  experiment,  the  ability  of  fluid  to  flow  through  cartilage  is  measured  for  a 
cartilage  specimen  in  a  fixed  position  (Mansour  et  al.  1976;  Torzilli  et  al.  1983;  Mak 
1986;  and  Gu  et  al.  1993).  The  stress  relaxation  experiment  involves  a  fast  cartilage  com¬ 
pression  (in  the  form  of  a  step  function).  The  stress  that  the  cartilage  specimen  imparts  on 
the  compressive  apparatus  is  then  measured  as  a  function  of  time  (Mow  et  al.  1980).  The 
shear  experiment  involves  fixing  a  mechanical  platen  or  indenter  to  the  cartilage  surface 
and  measuring  the  statics  and  dynamics  of  cartilage  by  imparting  a  torque  or  lateral  force 
(Parsons  and  Black  1977;  Spirt  et  al.  1989;  and  Zhu  et  al.  1993). 

The  present  work  focuses  on  developing  and  testing  a  useful  mechanical  model  of 
cartilage  compression.  The  most  commonly  accepted  existing  model,  developed  in  by 
V.C.  Mow  1980,  as  applied  to  the  uniaxial  compression  experiment  contains  a  set  of  cou¬ 
pled  partial  differential  equations  which  presents  a  challenge  for  many  cartilage  investi¬ 
gators  when  interpreting  experimental  results  in  light  of  contemporary  theory.  Mow’s 
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theoretical  model  may  be  appropriate  for  describing  cartilage  rheological  mechanisms 
from  a  continuum  mechanics  approach.  However,  from  a  practical  standpoint,  the  intui¬ 
tion  of  the  model  can  become  lost  in  its  complexity.  The  model  proposed  in  this  thesis  is 
comparatively  simple.  In  brief,  cartilage  is  treated  essentially  as  a  dialysis  bag  (i.e.  the 
mechanical  properties  of  cartilage  are  derived  primarily  from  osmotic  conditions  created 
by  the  cartilage  biochemical  composition.  The  fundamentals  of  the  model,  as  presented  in 
chapter  2,  are  relatively  easy  to  interpret.  An  additional  advantage  of  this  new  model  is 
that  it  easily  accommodates  anisotropy  and  inhomogeneity  as  addressed  further  in  appen¬ 
dix  B. 

An  interesting  aspect  to  the  model  is  that  water  and  sodium  ion  concentrations  are 
fundamental  to  cartilage  mechanical  properties.  This  fact  makes  magnetic  resonance  im¬ 
aging  (MRI)  of  cartilage  deformation  valuable  for  several  reasons.  MRI,  in  general,  is  an 
active  research  field  since  MRI  is  currently  the  most  promising  non-invasive  modality  to 
stage  osteoarthritis  progression.  Presented  here  is  MRI  conducted  in  a  multinuclear  fash¬ 
ion  to  noninvasively  map  cartilage  water  and  sodium  ion  content.  In  this  way,  MRI  pro¬ 
vides  an  excellent  means  for  testing  the  proposed  theoretical  model. 

1.4  Imaging  Cartilage  Compression 

As  mentioned  above,  magnetic  resonance  imaging  (MRI)  shows  promise  as  a 
non-invasive  modality  for  imaging  cartilage  changes,  since,  in  general,  MRI  is  able  to 
show  soft  tissue  contrast.  However,  current  attempts  at  using  MRI  to  measure  mild 
changes  in  articular  cartilage  have  encountered  inconclusive  results  (Hutton  et  al.  1995). 
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It  is  possible  that  an  understanding  of  the  MR  characteristics  of  articular  cartilage  during 
compression  will  lead  to  a  sensitive  non-invasive  indicator  of  disease  progression. 

The  earliest  of  the  cartilage  compression  studies  within  an  MRI  environment  was 
conducted  by  Lehner  et  al.  (1989).  This  experiment  employed  the  cylinder  and  plunger  of 
a  common  syringe  to  compress  a  glass  marble  against  the  surface  of  an  articular  cartilage 
specimen.  In  distinction  to  this  mechanical  approach,  Lusse  et  al.  (1995)  combined  MRI 
with  an  osmotic  compression  technique.  In  this  approach,  cartilage  specimens  were  en¬ 
closed  in  a  dialysis  membrane  then  placed  in  a  polyethylene  glycol  solution  of  high  os- 
molarity.  This  process,  in  effect,  dehydrated  the  cartilage  specimens.  Once  dehydrated, 
the  cartilage  specimens  were  removed  from  the  polyethylene  glycol  bath  then  placed  in 
an  MR  spectrometer  for  subsequent  measurement  of  MR  parameters. 

More  recently,  Rubenstein  et  al.  (1996)  performed  mechanical  compression  of  ar¬ 
ticular  cartilage  specimens  within  an  MRI  environment  using  a  compression  device  de¬ 
sign  based  on  an  earlier  description  at  the  1995  annual  meeting  of  the  International  Soci¬ 
ety  of  Magnetic  Resonance  in  Medicine  (Kaufman  et  al.  1995).  Further  experiments  ex¬ 
panded  upon  this  model  (Duvvuri  et  al.  1998,  Kaufman  et  al.  1999,  Reddy  et  al.  1999), 
several  of  which  constitute  the  basis  of  portions  of  the  work  presented  here. 

1.5  Thesis  Summary 

The  central  theme  of  this  thesis  involves  novel  methods  for  using  MRI  to  explore 
the  mechanical  properties  of  articular  cartilage  in  light  of  the  osteoarthritic  disease  proc¬ 
ess.  In  the  following  chapter,  MRI  is  used  as  a  tool  for  exploring  the  osmotic  ion  balance 

of  articular  cartilage,  which  imparts  the  majority  of  cartilage  resistance  to  mechanical 
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compression.  The  design  and  implementation  of  a  novel  MRI  compatible  cartilage  com¬ 
pression  device  is  then  described.  This  device  is  then  used  for  acquiring  MR  spectro¬ 
scopic  information  on  cartilage  specimens  during  mechanical  compression.  An  evaluation 
of  the  statistical  significance  of  compression  induced  spectroscopic  change  is  performed 
with  respect  to  an  in-vitro  osteoarthritic  model.  Spatial  mapping  of  the  spectroscopic  pa¬ 
rameters  is  then  performed  during  compression  to  investigate  the  spatial  redistribution  of 
MR  relaxation  characteristics  of  articular  cartilage  during  compression.  A  novel  method 
for  quantifying  the  dynamics  of  cartilage  compression  within  an  imaging  environment  is 
then  described  and  demonstrated.  The  results  of  this  method  are  then  used  to  infer  the 
compression-induced  dynamics  of  cartilage  fluid  permeability.  The  mechanical  predic¬ 
tions  of  the  osmotic  cartilage  model  are  formulated  and  confirmed  by  experiments.  An 
example  of  multidimensional  cartilage  deformation  is  then  demonstrated  and  described. 
Lastly,  a  technique  for  articular  cartilage  strain  mapping  during  compression  is  described. 

1.6  Thesis  Hypothesis  and  Strategy 

The  basic  working  hypothesis  for  the  present  work  is  that  ionic  equilibria  can  ac¬ 
curately  predict  cartilage  compressive  strength.  This  present  work  also  demonstrates  a 
new  technology:  the  ability  to  spatially  map  cartilage  mechanical  properties.  The  overall 
strategy  has  three  components.  First,  a  useful  deformation  theory  of  articular  cartilage  is 
developed.  Next,  a  device  is  developed  to  compress  cartilage  in  an  MRI  environment  and 
confirm  the  deformation  theory  with  experimental  data.  MRI  is  then  used,  in  a  novel  way, 
to  characterize  the  spatial  distribution  of  articular  cartilage  mechanical  properties. 
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Figure  1.1.  V.C.  Mow  model  for  cartilage  deformation  in  the  uniaxial  compression 
periment.  (From  Mow  et.  al.  1980.) 


Chapter 2 


A  Magnetic  Resonance  Imaging  Based  Demonstration 
of  Articular  Cartilage  Ionic  Equilibria 


2.1  Introduction 

As  mentioned  in  the  previous  chapter,  ionic  equilibria  is  an  important  component 
of  articular  cartilage  functionality.  Cartilage  is  a  stiff  tissue  with  strong  resistance  to  com¬ 
pression  while  being  primarily  composed  of  water.  Approximately  80%  of  cartilage  by 
weight  is  water.  The  stiffness  of  cartilage  results  from  the  water  being  contained  at  high 
pressure.  The  high  water  pressure  is  an  indirect  response  to  the  large  negative  fixed 
charge  density  associated  with  the  proteoglycan  component  of  the  cartilage  matrix.  It  has 
recently  become  popular  to  exploit  the  ionic  equilibria  of  cartilage  to  non-invasively 
probe  cartilage  proteoglycan  loss  via  charged  MRI  contrast  agents  (Bade  et  al.  1997; 
Bashir  et  al.  1997;  Bashir  et  al.  1996). 

2.2  Theory  of  cartilage  ionic  equilibria 

Consider  a  cartilage  specimen  equilibrated  in  an  infinitely  large  NaCl  bath.  To 
simplify  the  algebra,  all  concentrations  are  normalized  to  the  bathing  NaCl  concentration. 
The  cartilage  specimen  proteoglycan  content  imparts  an  internal  negative  fixed  charge 
density  in  the  cartilage  fluid  volume,  fi,  and  consequently  the  cartilage  fluid  volume  must 
have  a  sodium  ion  concentration  of  at  least/-  to  maintain  electrostatic  neutrality.  An  addi¬ 
tional  amount  of  sodium  content  will  equilibrate  into  the  cartilage  fluid  volume,  since 
NaCl  must  enter  the  cartilage  to  satisfy  the  chloride  ion  chemical  potential.  Thus, 
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n*  =  ft  +x, 


[2.1] 


where  n*  is  the  sodium  ion  concentration  within  the  cartilage  specimen  fluid  compart¬ 
ment,  is  the  fixed  charge  density  within  the  cartilage  specimen  fluid  compartment,  and 
x  is  the  amount  of  additional  NaCl  that  enters  the  cartilage  fluid  compartment.  A  general 
thermodynamic  condition  for  ionic  equilibria  between  cartilage  and  the  bath  is  that  the 
electrochemical  potential  of  every  component  be  equal  (Helfferich  1962),  i.e. 

—  =  4,  [2.2] 

n  c* 

where  n*  is  the  cartilage  fluid  sodium  ion  concentration,  n  is  the  bath  sodium  ion  concen¬ 
tration,  c*  is  the  cartilage  fluid  chloride  ion  concentration,  and  c  is  the  bathing  chloride 
ion  concentration.  Given  the  normalization  convention, 

n  =  c  =  1 ,  [2.3] 


combining  equations  2.1  and  2.3  yields 

n*c*  =  l.  [2.4] 


As  indicated  in  equation  2.1,  x  is  the  amount  of  NaCl  that  enters  the  cartilage 
from  the  bath  given  an  initial  condition  of  zero  cartilage  chloride  ion  content.  Thus, 

c*  =  x.  [2.5] 

Combining  equations  2.1, 2.4,  and  2.5  yields  the  quadratic  relationship 

(fi+x)x  =  1,  [2.6] 


which  has  the  solution 


x  = 


± 


[2.7] 
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Combining  equations  2.3, 2.5,  and  2.7  yields  the  chloride  ion  equilibrium  distribution,  d>, 
between  the  bath  and  the  cartilage  fluid  content 


[2.8] 


which  expressed  as  a  function  of  the  excluded  volume  fraction,  e,  and  the  fixed  charge 
density  within  the  total  cartilage  volume,/,  is 


<!>  =  _ 


/ 


•  +  Jl  + 


2{\-e)  ^  \2{\-e)) 


f 


[2.9] 


2.3  The  Gd(DTPA)  method  -  working  hypothesis 

Gadolinium  Diethylenetriaminepentaacetate,  Gd(DTPA),  has  captured  recent  in¬ 
terest  in  cartilage  research,  since  magnetic  resonance  imaging  (MRI)  techniques  can  po¬ 
tentially  map  its  concentration  to  high  spatial  resolution.  Gd(DTPA)  equilibration  pro¬ 
vides  one  avenue  of  investigation  toward  using  the  principles  of  ionic  equilibria  to  non- 
invasively  derive  the  proteoglycan  content  of  intact  articular  cartilage,  a  parameter  that 
may  be  indicative  of  osteoarthritic  progression  (Hardingham  et  al.  1990).  This  approach 
has  recently  been  established  in  an  in-vitro  setting  (Bashir  et  al.  1996),  and  extended  to 
in-vivo  investigations  (Bashir  et  al.  1997;  Trattnig  et  al.  1999).  To  use  this  method  one 
first  employs  standard  techniques  to  quantify  the  change  in  the  proton  spin-lattice  relaxa¬ 
tion  time  constant,  Ti,  associated  with  Gd(DTPA)  equilibration.  By  using  known  spin- 
lattice  relaxivities  of  Gd(DTPA),  the  Gd(DTPA)  equilibrium  between  the  cartilage  inte¬ 
rior  and  exterior  is  then  quantified.  Classical  Donnan  theory  of  compartmental  equilibria 
(Donnan  1924)  is  then  employed  to  calculate  the  negative  fixed  charge  density  that  would 
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result  in  the  measured  Gd(DTPA)  equilibrium.  The  cartilage  fixed  charge  density  is  then 
translated  into  cartilage  proteoglycan  content  by  using  the  commonly  accepted  anionic 
content  values  associated  with  common  cartilage  proteoglycan  aggregates  (Maroudas  et 
al.  1969). 

NMR  techniques  have  quantitatively  validated  the  Donnan  model  with  respect  to 
the  distribution  of  Na+  ions  in  cartilage  (Lesperance  et  al.  1992).  However,  several  factors 
can  cause  osmotic  systems  to  deviate  from  the  ideal  Donnan  model  (Helfferich  1962), 
and  previous  experiments  suggest  that  Gd(DTPA)  in  cartilage  may  not  obey  ideal  Don¬ 
nan  equilibria  when  the  cartilage  fixed  charge  density  is  greater  than  30  mM  (Donahue  et 
al.  1994).  The  objective  of  this  investigation  was  to  quantify  the  equilibria  of  Gd(DTPA) 
between  articular  cartilage  and  a  bathing  NaCl  solution,  and  thus  improve  our  under¬ 
standing  of  the  role  that  ionic  equilibria  has  on  cartilage  functionality  as  well  as  augment 
the  interpretation  of  Gd(DTPA)  related  investigations  of  cartilage  proteoglycan  content. 

2.4  Methods 

2.4.1  Overall  Protocol 

Cartilage  specimens  were  obtained  from  bovine  patellae.  All  specimens  were 
equilibrated  in  137  mM  phosphate  buffered  saline  (PBS).  After  equilibration,  one  dimen¬ 
sional  spin-lattice  relaxation  time  constant  (Ti)  maps  were  collected  along  an  arbitrary 
diameter  of  each  cartilage  disk.  Each  cartilage  specimen  was  then  bisected  along  a  diame¬ 
ter.  Half  of  each  specimen  was  equilibrated  in  137  mM  PBS  with  ImM  Gd(DTPA),  while 
the  other  half  of  each  specimen  was  equilibrated  in  varying  concentrations  of  PBS  with  1 

mM  Gd(DTPA).  The  specimen  halves  were  then  reunited  for  a  second  series  of  one  di- 
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mensional  Ti  maps,  along  the  diameter  perpendicular  to  specimen  bisection.  Two- 
dimensional  Ti  maps  were  also  acquired  to  evaluate  the  spatial  homogeneity  of  the 
Gd(DTPA)  associated  Ti  reduction. 

The  quantification  of  glycosaminoglycans  in  the  cartilage  was  conducted  spectro- 
photometrically  via  a  1,9-dimethylmethylene  blue  (DMMB)  assay  (Famdale  et  al.  1982). 

2.4.2  Cartilage  Specimen  Preparation 

Patellae  were  obtained  from  veal  calves  within  6  hours  of  sacrifice,  then  stored  on 
ice.  Cylindrical  cartilage  specimens,  8.0  mm  in  diameter  were  cut  in  the  same  flat  region 
of  each  patella.  Each  specimen  was  subsequently  removed  from  the  underlying  cartilage- 
bone  interface.  Upon  removal,  each  specimen  was  placed  in  137  mM  phosphate  buffered 
saline  (PBS)  (Sigma  P-4417)  at  0°C  with  pH  =  7.4.  The  mass  of  each  cartilage  specimen 
was  measured  periodically  after  placement  into  PBS  until  an  asymptotic  mass  value  was 
reached.  This  process  took  approximately  two  hours. 

2.4.3  Cartilage  Specimen  Equilibration 

All  cartilage  equilibrations  were  conducted  at  0°C.  Buffer  solutions  were  made  by 
making  four  serial  dilutions  to  an  8X  stock  of  PBS  (Sigma  P-4417).  The  resulting  solu¬ 
tions  were  8.0X,  4.0X,  2.0X,  1.0X,  and  0.5X  PBS.  The  contents  of  the  PBS  solution  at 
IX  were:  0.01  M  phosphate  buffer,  0.0027  M  potassium  chloride  0.137  M  sodium  chlo¬ 
ride,  with  a  pH  of  7.4.  Gd(DTPA)  (Berlex  NDC  50419-188-02  )  was  added,  to  each  of 
the  five  resulting  buffer  solutions,  to  1.0  mM.  After  collecting  imaging  data  on  the  carti¬ 
lage  specimens  without  Gd(DTPA),  the  specimens  were  placed  in  the  respective  50  mL 

solutions  that  contained  1 .0  mM  Gd(DTPA).  The  cartilage  specimens  were  equilibrated 
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in  these  solutions  for  a  period  of  420  min.  The  mass  of  each  cartilage  specimen  was 
measured  before  and  after  the  1.0  mM  Gd(DTPA)  equilibrations. 

2.4.4  Ti  Projections 

The  one-dimensional  Ti  (proton  spin-lattice  relaxation  time)  projection  maps 
were  acquired  in  a  2.0  Tesla  field.  A  local  gradient  set  was  used,  which  was  capable  of 
producing  8  Gauss/cm.  A  solenoidal  8  turn  1  cm  rf-coil  was  used.  The  direction  of  the 
projection  images  was  perpendicular  to  that  of  the  2.0  Tesla  field. 

A  series  of  projection  images  was  collected  to  generate  the  Ti  projection  maps. 
Each  projection  image  was  generated  via  a  90°-180°  spin  echo  sequence,  with  no  phase 
encoding  and  with  frequency  encoding  in  the  direction  of  the  projection.  To  map  Ti, 
steady  state  repetitive  projection  images  were  obtained  for  various  repetition  times.  For 
conditions  prior  to  Gd(DTPA)  incubation,  the  repetition  times  were  0.100  sec,  0.250  sec, 
0.500  sec,  1.000  sec,  1.500  sec,  2.000  sec,  4.000  sec,  and  6.000  sec.  For  conditions  after 
the  Gd(DTPA)  incubation,  the  repetition  times  were  0.070  sec,  0.100  sec,  0.250  sec, 
0.500  sec,  1.000  sec,  1.500  sec,  2.000  sec,  and  3.000  sec.  All  one  dimensional  spin-echo 
projection  images  employed  an  echo  time  of  20  msec. 

For  each  pixel,  the  image  intensity  from  the  one-dimensional  images,  at  steady- 
state  for  a  given  repetition  time,  was  fit  via  a  least  squares  technique  to  equation  2.10. 

I  =  A(l-e-kT«)+C  [2.10] 

where  /  is  the  steady-state  image  intensity,  k  is  1/T /,  TR  is  the  repetition  time,  A  is  the 
equilibrium  z-magnetization  corrected  for  T2  relaxation,  and  C  is  an  arbitrary  constant. 
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Prior  to  Gd(DTPA)  incubation,  the  cartilage  specimens  were  in  the  form  of  whole 
discs.  After  incubation  in  Gd(DTPA),  the  cartilage  specimens  were  in  the  form  of  half 
discs.  (See  the  ‘2.4.1  Overall  Protocol’  section.)  The  one-dimensional  projection  maps 
prior  to  Gd(DTPA)  incubation  were  taken  along  an  arbitrary  diameter.  After  Gd(DTPA) 
incubation,  all  one-dimensional  projections  were  taken  with  the  two  halves  of  each 
specimen  placed  adjacent  to  each  other  along  their  diameter  of  bisection,  and  oriented  in 
the  gradient  set  such  that  each  one-dimensional  projection  image  was  collected  along  a 
direction  perpendicular  to  the  diameter  of  bisection.  All  cartilage  specimens  were  ori¬ 
ented  such  that  the  articular  surface  was  parallel  to  the  direction  of  the  2.0  Tesla  magnetic 
field. 

2.4.5  Two-Dimensional  Ti  Maps 

Two-dimensional  slice  selective  Tj  maps  were  collected  to  detect  the  possibility 
of  an  uneven  Gd(DTPA)  distribution.  Each  slice  was  collected  along  the  same  diameter 
as  that  of  the  one-dimensional  projection  images.  Slice  thickness  was  2.0  mm.  The  two 
dimensional  maps  were  generated  from  a  series  of  two  dimensional  90°-180°  spin-echo 
images.  Each  two-dimensional  image  employed  a  20  msec  echo  time.  Frequency  encod¬ 
ing  was  in  the  direction  perpendicular  to  the  cartilage  articular  surface.  Phase  encoding 
was  in  the  direction  parallel  to  the  cartilage  articular  surface.  Different  repetition  times 
were  employed  to  create  the  individual  images  with  which  the  Ti  map  was  generated. 
The  repetition  times  used  were  100  msec,  250  msec,  400  msec,  750  msec,  1500  msec, 
and  3000  msec.  Each  pixel  of  the  two-dimensional  images  was  fit  to  equation  2.10  via  a 
least  squares  technique. 
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2.5  Results 


2.5.1  Tj  Projections  of  Cartilage  Specimens 

The  Ti  measurements  of  the  cartilage  specimens  are  shown  in  table  2.1.  For  each  carti¬ 
lage  specimen,  the  Tj  maps  were  homogeneous  with  random  spatial  fluctuations  to  the 
extent  indicated  by  the  standard  deviations  of  the  Ti  values.  The  Ti  spatial  fluctuations 
were  confirmed  to  be  random  by  demonstrating  that  a  given  Ti  spatial  pattern  was  not 
reproducible  when  repeating  the  Ti  mapping  procedure  on  a  cartilage  specimen  several 
times.  The  histograms  of  Ti  values  within  each  cartilage  specimen  were  approximately 
gaussian,  as  shown  by  a  representative  Ti  histogram  in  figure  2.1.  Each  Ti  projection 
contained  approximately  230  pixels.  The  homogeneity  of  the  Ti  maps,  in  a  direction  par¬ 
allel  to  the  articular  surface,  is  demonstrated  by  the  lack  of  non-random  Ti  variations  in  a 
typical  one-dimensional  projection  map  shown  in  figure  2.2. 

2.5.2  Gd(DTPA)  Induced  Tj  Depression 

The  Ti  measurements  indicated  that  the  Gd(DTPA)  associated  Ti  reduction  in  the  carti¬ 
lage  specimens  was  more  pronounced  when  the  specimens  were  equilibrated  in  higher 
PBS  concentrations.  The  means  and  standard  deviations  of  the  Ti  values,  calculated  from 
the  Ti  projection  maps,  for  each  specimen  half  are  shown  in  table  2.2.  Gd(DTPA)  con¬ 
centrations  in  the  cartilage  specimen  halves  were  calculated  via  equation  2.1 1 

C211] 

where  g  is  the  total  cartilage  Gd(DTPA)  concentration  in  mM,  R  is  the  relaxivity  of 
Gd(DTPA)  in  (mM-sec)'1,  rt*  is  the  cartilage  proton  spin-lattice  relaxation  rate  after 
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Gd(DTPA)  equilibration,  and  n  is  the  cartilage  proton  spin-lattice  relaxation  rate  prior  to 
Gd(DTPA)  equilibration.  We  employed  a  previously  reported  estimate  for  Gd(DTPA) 
relaxivity,  4.0  (mM-sec)'1(Donahue  et  al.  1994).  The  results  of  these  calculations  are 
shown  in  figure  2.3. 

Figure  2.4  illustrates  the  asymptotic  limit  of  Gd(DTPA)  equilibria  as  a  function  of 
the  ratio  of  the  bathing  NaCl  concentration  and  the  cartilage  fixed  charge  density.  The 
experimental  data  has  a  reasonable  fit  (r  =  0.99516)  to  a  Donnan  model  that  incorporates 
a  cartilage  volume  excluded  to  Gd(DTPA)  penetration.  This  model  is  described  by  equa¬ 
tions  2.12  and  2.13, 

^  =  (l-*)bz  [2.12] 
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where  d>  is  the  cartilage  chloride  ion  equilibrium  distribution,  g*  is  the  Gd(DTPA)  con¬ 
centration  inside  the  total  cartilage  volume,  g  is  the  Gd(DTPA)  concentration  outside  the 
cartilage, /is  the  total  cartilage  negative  fixed  charge  density  as  a  fraction  of  the  bathing 
NaCl  concentration,  e  is  the  fraction  of  cartilage  volume  excluded  from  Gd(DTPA)  pene¬ 
tration,  and  Z  is  the  apparent  charge  of  Gd(DTPA).  Equation  2.12  is  an  analog  of  equa¬ 
tion  2.2,  and  equation  2.13  is  derived  in  section  2.2.  Least  squares  fit  of  the  model  to  the 
experimental  results  are  shown  in  the  figure  for  Z  =  -2  and  Z  =  -1.  An  apparent 
Gd(DTPA)  charge  of  -1  provides  a  much  better  fit  than  the  does  an  apparent  charge  of -2. 
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When  the  apparent  Gd(DTPA)  charge  was  used  as  a  fit  parameter,  the  best  fit  yielded  an 
excluded  volume  of  20%  and  Gd(DTPA)  effective  charge  of -0.92. 

2.5.3  Two-dimensional  Ti  Mapping 

The  one  dimensional  maps  of  figure  2.3  show  a  uniform  Ti  distribution  within 
each  specimen  half,  which  indicates  a  likelihood  that  there  was  enough  time  for 
Gd(DTPA)  to  fully  equilibrate.  However,  a  uniform  two  dimensional  Ti  map  lends  fur¬ 
ther  support  to  this  notion,  and  reduces  the  possibility  that  the  Gd(DTPA)  diffusion  was 
incomplete  in  a  direction  perpendicular  to  the  cartilage  specimen’s  surface.  Figure  2.6 
shows  a  slice  selective  two  dimensional  Ti  map  generated  from  multiple  two  dimensional 
spin  echo  images  of  a  cartilage  specimen  after  bisection  and  Gd(DTPA)  equilibration. 
Each  half  of  the  cartilage  specimen  has  a  uniform  Ti  distribution  with  a  similar  level  of 
Ti  noise  as  was  seen  in  the  one  dimensional  Ti  projection  images  of  figure  2.3. 

2.6  Discussion 

Several  MRI  techniques  to  non-invasively  measure  proteoglycan  content  have 
been  explored,  such  as  the  use  of  sodium  NMR  spectroscopy  (Lesperance  et  al.  1992),  in- 
vivo  sodium  imaging  of  articular  cartilage  (Reddy  et  al.  1998;  Shapiro  et  al.  2000),  and 
the  use  of  contrast  agents  (Bacic  et  al.  1997),  including  Gd(DTPA)  (Bashir  et  al.  1996; 
Bashir  et  al.  1997).  The  present  study  has  used  MRI  to  quantify  the  equilibrium  of 
Gd(DTPA)  between  articular  cartilage  and  a  bathing  saline  solution,  and  tested  the  hy¬ 
pothesis  that  the  Gd(DTPA)  distribution  is  governed  by  ideal  Donnan  equilibria.  A  quan¬ 
titative  understanding  of  Gd(DTPA)  equilibria  in  cartilage  is  essential  to  a  meaningful 

interpretation  of  associated  derivations  of  cartilage  proteoglycan  content. 
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The  expectation  that  the  fixed  charge  density  associated  with  cartilage  proteogly¬ 
can  content  reduces  Gd(DTPA)  transport  into  the  cartilage  interior  was  experimentally 
realized.  The  results  were  consistent  with  a  Donnan  model  that  included  a  provision  for  a 
cartilage  volume  inaccessible  to  Gd(DTPA)  entry.  The  modified  Donnan  model  derived 
here  is  different  than  the  commonly  accepted  model  (Bashir  et  al.  1996)  in  two  respects. 
First  the  present  model  has  an  explicit  incorporation  of  a  cartilage  volume  inaccessible  to 
Gd(DTPA)  equilibration.  Second,  the  present  model  does  not  dictate  that  the  effective 
charge  of  the  Gd(DTPA)  complex  is  -2. 

The  prediction  that  the  Gd(DTPA)  charge  is  considerably  less  negative  than  -2 
under  physiological  conditions  is  not  inconsistent  with  previous  equilibrium  investiga¬ 
tions  of  metal  chelating  agents.  DTPA  alone  has  five  potentially  ionizable  groups  with  a 
wide  range  of  pK  values:  1.79,  2.56,  4.42,  8.76,  and  10.42  (Werner  and  Grunder  1999), 
and  it  is  common  for  DTPA  metal  chelates  to  be  protonated  at  neutral  pH  (Durham  et  al. 
1958).  It  is  not  clear  what  the  protonation  equilibrium  properties  of  Gd(DTPA)  are  in  the 
cartilage  interior,  and  thus  further  experimentation  would  be  required  to  confirm  the  ef¬ 
fective  Gd(DTPA)  charge  in  cartilage.  The  effective  charge  estimated  in  figure  5  should 
be  considered  a  purely  empirical  result. 

The  fit  to  the  Donnan  model  indicates  that  the  cartilage  volume  fraction  available 
to  Gd(DTPA)  is  approximately  80%.  The  major  cartilage  biopolymers,  collagen  and 
agrecan,  constitute  a  substantial  cartilage  volume  fraction:  approximately  20%  as  meas¬ 
ured  in  our  laboratory.  Thus  the  accessible  cartilage  volume  fraction  unavailable  to 
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Gd(DTPA)  penetration  is  likely  to  be  coincident  to  the  volume  occupied  by  the  major 
cartilage  biopolymers. 

The  results  of  this  investigation  indicate  that  efforts  to  quantify  cartilage  proteoglycan 
content  from  Gd(DTPA)  equilibria  should  take  into  account  the  apparent  charge  of  the 
Gd(DTPA)  complex  and  the  volume  effects  of  the  cartilage  solid  content.  Cartilage  solid 
content  is  an  especially  important  parameter  in  this  consideration  since  cartilage  solid 
content  may  vary  during  pathological  conditions. 
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Table  2.1.  Ti  measurements  (at  2  Tesla)  of  the  five  cartilage  specimens  equilibrated  in 
137  mM  PBS  without  Gd(DTPA).  Standard  deviations  were  obtained  by  considering 
each  pixel  of  the  Ti  projection  maps  to  be  one  element  of  a  population  of  Ti  measure¬ 
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Figure  2.1.  Distribution  of  Ti  values  of  a  representative  cartilage  specimen  equilibrated 
in  137  mM  PBS  without  Gd(DTPA). 
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Figure  2.2.  One-dimensional  Ti  projection  map  of  a  representative  cartilage  specimen 
equilibrated  in  137  mM  PBS  without  Gd(DTPA).  The  thickness  of  one  pixel  is  approxi¬ 
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Table  2.2.  Ti  measurements  of  the  five  cartilage  specimens  equilibrated  in  PBS  solutions 
of  various  NaCl  concentrations  with  1.0  mM  Gd(DTPA).  Standard  deviations  were  ob¬ 
tained  by  considering  each  pixel  of  the  Ti  projection  maps  to  be  one  element  of  a  popula¬ 
tion  of  Ti  measurements. 
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[FCD]/[NaCI] 

Figure  2.3.  Partition  coefficient  of  Gd(DTPA)  transport  between  the  cartilage  interior 
and  the  bathing  medium  as  a  function  of  fixed  charge  density  normalized  to  the  NaCl 
concentration  in  the  bathing  medium. 
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NaCI/FCD 

Figure  2.4.  Gd(DTPA)  transport  partition  coefficient  as  defined  in  figure  2.3  plotted  as  a 
function  of  the  multiplicative  inverse  of  the  fixed  charge  density.  The  solid  lines  are  least 
squares  fits  of  equation  2.14  to  the  experimental  results.  The  fits  were  performed  for  two 
apparent  charges  (Z  =  -1  and  Z  =  -2)  as  indicated.  The  least  squares  fit  when  both  the 
charge  and  excluded  volume  were  considered  parameters  yielded  an  excluded  volume  of 
20%  and  an  effective  Gd(DTPA)  charge  of  -0.92. 


27 


Figure  2.5.  Slice  selective  two-dimensional  Ti  maps  of  a  cartilage  specimen.  The  right 
and  left  portions  of  the  figure  correspond  to  halves  of  a  cartilage  disk  after  bisection. 
Each  half  of  the  cartilage  disc  was  equilibrated  in  1.0  mM  Gd(DTPA)  and  saline  concen¬ 
trations  as  indicated  in  the  figure. 
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Chapter  3 


A  Method  for  Compressing  Articular  Cartilage 
Within  a  Magnetic  Resonance  Imaging  Environment 

3.1  Introduction 

The  normal  function  of  articular  cartilage  involves  compression.  As  articular  car¬ 
tilage  is  compressed,  the  fluid  contained  within  the  cartilage  matrix  is  extruded  into  the 
joint  space.  The  extruded  fluid  lubricates  the  contact  points  of  opposing  bone  interfaces 
(McCutchin  1972).  It  is  possible  that  cartilage  that  is  too  soft  will  compress  to  an  unnatu¬ 
ral  extent,  and  provide  a  mechanical  advantage  for  natural  shear  forces  to  disrupt  the  col¬ 
lagen  framework.  Additionally,  soft  cartilage  may  not  be  able  to  recover  its  original 
shape  quickly  enough  to  continue  providing  fluid  to  a  bone  interface  during  repeated 
loading.  To  test  hypotheses  about  cartilage  physiology  and  pathophysiology,  it  is  useful 
to  develop  non-invasive  imaging  techniques  for  observing  the  functional  performance  of 
articular  cartilage.  We  have  developed  a  technology  which  helps  explore  the  magnetic 
resonance  imaging  properties  of  articular  cartilage  specimens  while  simultaneously  ac¬ 
quiring  information  regarding  deformation.  The  design  of  a  magnetic  resonance  imaging 
compatible  cartilage  compression  device  is  presented. 

3.2  Design 

The  pressure  cell,  shown  in  figure  3.1,  consists  of  a  pneumatically  activated  ce¬ 
ramic  piston  confined  in  a  cylindrical  chamber.  The  cartilage  was  positioned  in  the  pres- 
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sure  cell  with  its  articular  surface  in  contact  with  a  porous  glass  filter  (Coming  31001- 
10F).  The  filter  led  to  a  drain,  which  allowed  extruded  cartilage  fluid  to  leave  the  cell. 

The  cylinder  and  cap  of  the  pressure  cell  were  made  of  Delrin  (Dupont).  This 
polymer  exhibits  no  proton  NMR  signal,  is  nonmetallic,  has  an  electrical  resistivity  on 
the  order  of  1014£2-cm  and  a  dissipation  factor  of  0.0048  at  106  Hz.  The  piston  was  made 
of  a  machinable  ceramic,  Macor  (Coming).  This  ceramic  has  an  electrical  resistivity  on 
the  order  of  1 016  Q-cm  and  a  loss  tangent  of  10 ' 3  at  1  kHz  and  104  at  8.5  GHz.  The 
pore  size  of  the  glass  filter  was  chosen  such  that  it  was  impermeable  to  the  solid  matrix  of 
the  cartilage  but  freely  permeable  to  the  cartilage  fluid.  The  compressed  air  that  pneu¬ 
matically  activated  the  piston  was  prevented  from  traveling  around  the  sides  of  the  piston 
by  a  Teflon  seal.  Pressure  was  applied  via  a  compressed  air  tank,  and  was  measured  by  a 
digital  pressure  gauge  (Omega  DPG-500)  on  the  regulator.  A  solenoid  radio-frequency 
coil  was  wrapped  around  the  cylinder  of  the  pressure  cell.  The  pressure  cell  was  placed 
inside  an  1 1  cm  diameter  local  gradient  set  capable  of  generating  8G/cm.  The  cell  and 
gradient  set  were  then  placed  in  the  1  meter  bore  of  a  2  Tesla  Oxford  whole-body  imag¬ 
ing  magnet  interfaced  to  a  custom-built  NMR  spectrometer.  The  functionality  of  the 
compression  cell  was  tested  by  performing  several  experiments  on  cartilage  plugs  as  de¬ 
scribed  below. 

3.3  Experimental  protocol 

Cartilage  specimens  were  obtained  from  the  patella  of  bovine  calves  (ages  ranged 

from  6  months  to  a  year)  within  10  hours  from  sacrifice.  The  cartilage  specimens  were 

stored  and  treated  slightly  differently  then  were  the  specimens  used  in  chapter  2.  Each 
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specimen  was  stored  at  -20°C  for  a  period  no  longer  than  20  days.  Cartilage  specimens 
were  harvested  from  the  flat  portion  of  the  articular  surface  of  the  patella.  With  a  small 
spatula,  each  disc  was  peeled  from  the  cartilage-bone  interface.  Prior  to  insertion  into  the 
pressure  cell,  each  specimen  was  equilibrated  in  10  mL  of  phosphate  buffered  saline 
(Sigma  P-4417)  at  room  temperature  for  30  minutes. 

In  an  effort  to  simulate  articular  cartilage  in  an  early  arthritic  state,  we  exposed 
several  cartilage  specimens  to  a  moderate  proteolytic  regime.  Since  it  has  been  hypothe¬ 
sized  that  early  stages  of  osteoarthritis  involve  the  loss  of  proteoglycan  content  selec¬ 
tively,  we  chose  a  proteolytic  agent  that  would  selectively  deplete  the  proteoglycan  con¬ 
tent  of  the  cartilage  while  leaving  the  collagen  network  relatively  intact.  The  proteolytic 
enzyme,  trypsin,  is  known  to  cleave  the  core  proteins  of  the  proteoglycans  while  having 
little  effect  on  the  type-II  collagen  fibers  which  comprise  the  collagenous  network  of  ar¬ 
ticular  cartilage.  The  digested  cartilage  specimens  referred  to  in  this  chapter  were  gener¬ 
ated  by  incubating  each  cartilage  specimen  in  a  10  mL  volume  of  1  mg/mL  trypsin- 
(Sigma  T-8253)/phosphate-buffered-saline-(Sigma  P-4417)  solution  for  75  minutes  at 

37°C  with  slight  agitation. 

Slice  selective  images  of  each  cartilage  specimen  in  a  plane  coincident  with  a  di¬ 
ameter  of  the  cartilage  disk  and  perpendicular  to  the  cartilage  surface  were  acquired. 
Slice  thickness  was  2  mm.  In  each  case  the  cartilage  specimen  was  oriented  such  that  its 
articular  surface  was  parallel  to  Bq,  the  main  magnetic  field.  All  two-dimensional  images 

were  acquired  with  a  standard  90°-180°  spin  echo  rf-pulse  sequence.  A  10  kHz  band¬ 
width  was  used.  The  orientations  of  frequency  and  phase  encoding  are  indicated  in  figure 
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3.1.  Each  image  employed  a  256  x  256  array  of  pixels.  All  images  were  acquired  using  a 
2  Tesla  full  body  imaging  magnet  integrated  with  a  custom  built  spectrometer.  Local 
gradients  that  enclose  the  pressure  cell  were  constructed  that  provide  up  to  8  gauss  per 
cm.  All  rf  pulses  were  applied  via  a  solenoidal  coil  tuned  to  the  proton  frequency  having 
6  turns  and  a  diameter  of  2.5  cm.  Imaging  data  were  analyzed  using  the  Interactive  Dis¬ 
play  Language  interpreter  on  a  Silicon  Graphics  workstation. 

Pressure  was  applied  to  the  cartilage  in  the  form  of  a  step  function.  Prior  to  pres¬ 
sure  application,  a  series  of  two-dimensional  90°-180°  spin-echo  images  were  acquired. 
The  compression  was  driven  by  a  pressure  of  0.690  MPa  (lOOpsi).  After  the  80  minutes 
of  constant  pressure  application,  the  pressure  of  0.690  MPa  was  maintained  for  a  subse¬ 
quent  set  of  image  acquisitions.  The  cartilage  fluid  which  emanated  from  the  pressure  cell 
drain  was  assayed  for  the  presence  of  proteoglycan  content  via  a  dimethylmethylene  blue 
spectrophotometric  assay  (Famdale  et  al.  1982). 

3.4  Results 

Figures  3.2  and  3.3  illustrate  compression  of  both  normal  and  proteoglycan  depleted  car¬ 
tilage  specimens.  The  signal  to  noise  ratio  of  these  images  is  approximately  7.  On  aver¬ 
age,  the  normal  cartilage  specimens  compressed  by  approximately  38%.  The  proteogly¬ 
can  depleted  cartilage  specimens  compressed  by  approximately  50%.  The  result  that  the 
proteoglycan-depleted  specimens  were  less  resistant  to  mechanical  compression  is 
consistent  with  the  notion  that  a  reduction  in  cartilage  fixed  charge  density  leads  to  a 
reduced  osmotic  pressure  within  the  cartilage  interior.  Although  the  viscoelastic  behavior 
of  articular  cartilage  is  not  simple,  a  description  of  the  equilibrium  resistance  of  articular 
cartilage  as  a  compressed  ideal  spring  proves  a  first  approximation  of  cartilage  com- 


lage  as  a  compressed  ideal  spring  provides  a  first  approximation  of  cartilage  compressive 
modulus  as  defined  by  equation  3.1 

G  —  Ee  [3.1] 

where  a  is  compressive  pressure,  E  is  compressive  modulus,  and  e  is  compressive  strain. 
In  this  analysis,  the  compressive  modulus  of  normal  articular  cartilage  is  approximately 
1.8  MPa,  and  the  compressive  modulus  of  proteoglycan  depleted  articular  cartilage  is  ap¬ 
proximately  1.38  MPa.  The  proteoglycan  depletion,  in  this  case,  caused  a  23%  loss  of 
compressive  resistance.  The  compressive  moduli  are  presented  here  with  the  caveat  that, 
according  to  prior  literature  (Mow  et  al.,  1980),  these  compressive  events  most  likely  ex¬ 
ceeded  the  compressive  elastic  limit  of  articular  cartilage,  and  thus  the  values  of  com¬ 
pressive  moduli  are  probably  overestimates  of  what  would  be  found  for  milder  compres¬ 
sive  events. 

It  should  be  noted  that  it  is  difficult  to  interpret  results  regarding  the  cartilage  ar¬ 
ticular  surface,  since  substantial  image  artifacts  exist  at  the  boundary  between  the  carti¬ 
lage  specimen  and  the  porous  glass  filter.  The  artifacts  present  as  localized  regions  of 
bright  signal  intensity  and  a  non-smooth  articular  surface. 

The  results  of  the  proteoglycan  assays  of  extruded  cartilage  fluid  indicated  no 
measurable  proteoglycan  content,  indicating  that  cartilage  compressive  events  involve 
minimal  extrusion  of  solid  content. 

3.5  Discussion 

The  functionality  of  a  novel  cartilage  compression  device  has  been  demonstrated. 

The  effective  compression  of  small  articular  cartilage  specimens  was  shown  for  both 
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normal  cartilage  specimens  and  cartilage  specimens  softened  by  selective  enzymatic  deg¬ 
radation.  The  degree  of  cartilage  compression  as  a  response  to  a  pressure  of  0.69  MPa  is 
consistent  with  previous  biomechanical  results  reported  in  the  literature  (Mow  et.  al. 
1980).  The  major  challenge  with  the  compression  device  is  that  small  yet  significant  im¬ 
aging  artifacts  are  present  at  the  cartilage  to  glass  filter  interface.  The  size  of  the  artifacts 
is  several  orders  of  magnitude  larger  than  the  roughness  of  the  filter  surface.  One  likely 
source  for  the  imaging  artifacts  is  the  possibility  of  magnetic  susceptibility  effects  result¬ 
ing  from  trapped  air  in  the  filter.  During  compression,  the  presence  of  image  artifacts 
substantially  reduce. 
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Figure  3.1.  Cross  sectional  diagram  of  MR  compatible  pressure  cell.  Compressed 
air  exerts  a  pressure  on  a  piston  that  is  confined  in  a  cylindrical  shaft.  The  piston, 
when  responding  to  the  pressure  from  the  compressed  air,  pushes  on  the  specimen 
that  is  positioned  between  the  cylinder  and  a  Pyrex  filter,  which  allows  cartilage 
fluid  to  extrude  from  the  cell. 
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Figure  3.2.  Two-dimensional  slice  selective  90°-180°  spin  echo  images  of  a  normal  ar¬ 
ticular  cartilage  specimen  prior  to  and  during  compression.  Cartilage  orientation  is  such 
that  the  articular  surface  is  on  the  bottom.  The  cartilage  specimen  is  compressed  from 
above  by  a  ceramic  piston  onto  a  porous  glass  filter.  Imaging  artifacts  from  the  cartilage 
surface  to  glass  filter  interface  are  visible  in  both  images,  and  are  more  apparent  in  the 
image  prior  to  compression. 


PG  Depleted  Specimen 


Figure  3.3.  Two  dimensional  slice  selective  90°-180°  spin  echo  images  of  a  proteoglycan 
depleted  articular  cartilage  specimen  prior  to  and  during  compression.  Cartilage  orienta¬ 
tion  is  such  that  the  articular  surface  is  on  the  bottom.  The  cartilage  specimen  is  com¬ 
pressed  from  above  by  a  ceramic  piston  onto  a  porous  glass  filter.  Imaging  artifacts  from 
the  cartilage  surface  to  glass  filter  interface  are  visible  in  both  images,  and  are  more  ap¬ 
parent  in  the  image  prior  to  compression. 
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Chapter  4 


Spectroscopic  Magnetic  Resonance  Imaging  Parameters 
of  Articular  Cartilage  During  Compression 

4.1  Purpose  and  objective 

Chapter  1  introduced  several  investigations  that  have  used  proton  MRI  to  observe 
cartilage  image  dynamics  during  compression  (Lehner  et  al.  1989;  Rubenstein  et  al. 
1996).  These  studies  could  only  report  qualitative  effects  because  MR  relaxation  time 
constants  were  not  measured.  In  the  present  experiment  we  use  standard  inversion  recov¬ 
ery  and  the  standard  (90°-180°)  spin-echo  technique  to  obtain  proton  MR  relaxation  char¬ 
acteristics  of  cartilage  during  mechanical  compression.  We  also  examined  the  effects  of 
enzymatic  proteoglycan  depletion  on  the  observed  compression  induced  trends  in  relaxa¬ 
tion  times. 

In  addition  to  spectroscopic  measurements,  the  MR  relaxation  times  of  articular 
cartilage  during  compression  are  studied  in  more  detail  by  constructing  two-dimensional 
maps  of  cartilage  MR  proton  relaxation  time  constants.  The  two-dimensional  maps  pre¬ 
vent  the  measurement  of  cartilage  MR  relaxation  characteristics  from  being  influenced  by 
water  signal  emanating  from  within  the  porous  Pyrex  filter  of  the  MR  compatible  com¬ 
pression  device.  Furthermore,  the  MR  relaxation  time  constant  maps  permit  the  observa¬ 
tion  of  spatial  variations  in  cartilage  MR  relaxation  properties. 

4.2  Experimental  protocol 

4.2.1  Cartilage  Specimens 
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Cylindrical  cartilage  plugs  (1  cm  in  diameter)  were  harvested,  each  from  the  flat 
surface  of  a  bovine  patella.  All  patellae  were  stored  at  -20°C  prior  to  the  experiment. 
Each  cartilage  specimen  was  sectioned  with  a  diamond  core  drill  then  peeled  from  the 
cartilage-bone  interface  with  a  spatula.  The  specimens  were  subsequently  inspected  to 
ensure  that  there  was  no  visible  damage  and  that  all  cartilage  layers  remained  intact.  Each 
specimen  was  equilibrated  in  phosphate  buffered  saline  (pH  =  7.4)  for  30  minutes. 

Several  cartilage  specimens  were  exposed  to  enzymatic  digestion.  Since  it  has 
been  hypothesized  that  early  stages  of  osteoarthritis  involve  the  selective  loss  of  pro¬ 
teoglycan  content,  we  chose  a  proteolytic  agent  that  would  deplete  the  proteoglycan  con¬ 
tent  of  the  cartilage  while  leaving  the  collagen  network  relatively  intact.  The  digested  car¬ 
tilage  specimens  were  generated  by  incubating  each  cartilage  specimen  in  5ml  of  1 .0  mg 
trypsin  /ml  phosphate  -buffered  saline  solution  for  60  minutes.  The  end  stage  digestion 
media  was  analyzed  for  proteoglycan  content  via  a  spectrophotometric  assay  (Famdale  et 
al.  1982).  In  each  case,  the  proteoglycan  concentration  in  the  media  was  approximately 
1.3  mg/ml  which  corresponds  to  an  approximate  30%  proteoglycan  depletion. 

4.2.2.  Cartilage  Compression 

The  cartilage  compression  device  and  rf-coil  design  were  the  same  as  described  in 
chapter  3.  Each  cartilage  compression  event  consisted  of  a  constant  pressure  of  0.69  MPa 
was  employed  for  a  duration  of  180  minutes. 

4.2.3.  NMR  Relaxation  Time  Constant  Measurements 

The  standard  inversion  recovery  technique  was  employed  to  obtain  proton  Ti 

measurements.  For  normal  cartilage  specimens,  the  inversion  time  was  varied  by  step 
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times  of  200  msec  before  compression  and  100  msec  during  compression.  For  the  pro¬ 
teoglycan  depleted  cartilage  specimens,  the  inversion  time  varied  by  step  times  of  250 
msec  before  compression  and  125  msec  during  compression.  In  each  case  30  free  induc¬ 
tion  decays  (fids)  were  acquired.  Each  fid  contained  1024  data  points  at  a  sampling  rate  of 
0.1  msec. 

To  measure  proton  T2,  the  standard  Hahn  spin-echo  technique  (Hahn  1950)  was 
employed,  as  modified  by  Carr  and  Purcell.  See  “Method  A”  within  (Carr  and  Purcell, 
1954).  For  normal  cartilage,  the  echo  time  ( tau )  was  varied  by  25  msec  before  compres¬ 
sion  and  10  msec  during  compression.  For  proteoglycan  depleted  cartilage,  tau  was  var¬ 
ied  by  20  msec  before  compression  and  8  msec  during  compression.  Spectra  at  30  time 
points  were  acquired  for  each  T2  measurement.  For  each  time  point,  200  fids  were  aver¬ 
aged  to  obtain  a  spectrum. 

4.2.4.  MR  Relaxation  Time  Constant  Mapping 

We  constructed  maps  of  the  spin-lattice,  Ti,  and  spin-spin,  T2,  relaxation  times  for 
each  position  of  the  cartilage  imaging  slice  corresponding  to  an  image  pixel.  In  general, 
the  relaxation  time  maps  were  constructed  from  fitting  equations  4.1  and  4.2 


I  =  A(l-exp(-TRJTl))  +  B 

[4.1] 

I  =  Cexp(-TE/T2) 

[4.2] 

where  I  is  image  intensity,  TR  is  the  repetition  time,  TE  is  the  echo  time,  and  A,B,C  are 
arbitrary  constants,  to  a  series  of  90°-180°  spin-echo  images  employing  different  echo 
times  and  different  repetition  times.  Five  images  were  used  to  construct  each  Ti  map.  All 
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were  acquired  with  an  echo  time  of  20  msec,  and  had  the  following  repetition  times:  0.5 
sec,  1.0  sec,  1.5  sec,  2.0  sec,  and  2.5  sec.  Six  images  were  used  to  construct  each  T2  map. 
All  were  acquired  with  a  repetition  time  of  2.0  sec,  and  had  the  following  echo  times:  20 
msec,  30  msec,  40  msec,  50  msec,  60  msec,  and  70  msec.  The  collected  images  were  fit 
to  equations  4.1  and  4.2  using  a  gradient  expansion  algorithm  for  a  non-linear  least 
squares  fit.  It  is  important  to  note  that  this  technique  for  mapping  the  T2  relaxation  time 
constant  is  not  independent  of  diffusion  effects,  and  thus  the  T2  values  that  we  report  de¬ 
viate  from  true  T2  values  since  our  method  has  not  incorporated  a  diffusion- weighting 
factor  (Neeman  et  al.  1990;  Brandi  et  al.1994).  As  early  as  1954,  Carr  and  Purcell  have 
demonstrated  the  contribution  of  diffusion  when  observing  T2  decay  with  a  series  of  indi¬ 
vidual  spin  echo  events  to  be 

Myd  =  -M0exp[fDt3/12n2],  [4.3] 

where  Myd  is  the  contribution  of  diffusion,  y  is  the  gyromagnetic  ratio,  D  is  the  diffusion- 
constant,  t  is  duration  that  a  spin  ensemble  experiences  a  gradient,  and  n  is  the  number  of 
echoes  (Carr  et  al.  1954).  Considering  our  experimental  parameters  of  gradients  no  larger 
than  8  gauss/cm,  and  using  a  diffusion  constant  of  10’5  cm2/sec,  we  estimate  that  diffu¬ 
sion  can  influence  our  T2  measurements  by  a  factor  of  approximately  10%. 

In  addition,  given  the  fixed  set  of  repetition  times,  our  method  for  mapping  Ti  has 
an  inherent  limitation  on  the  range  of  Ti  values.  This  limitation  is  especially  important 
when  considering  Ti  values  near  the  lowest  repetition  time  of  0.5  sec.  A  calibration  of 
our  Ti  mapping  technique  with  copper  sulfate  doped  agarose  phantoms  of  known  Ti  and 
cartilage  specimens  (where  we  measured  bulk  Ti  with  an  inversion  recovery  technique) 


41 


indicated  that  our  Ti  mapping  procedure  yielded  Ti  values  accurate  to  within  an  ap¬ 
proximate  50%  range. 

4.3  Results 

4.3.1  Spectrocopic  Measurements 

Figures  4. 1-4.3  show  typical  curve  fits  when  calculating  relaxation  time  constants. 
Table  4.1  presents  a  summary  of  the  data  obtained  regarding  the  effects  of  compression 
and  proteoglycan  depletion  on  cartilage  proton  and  sodium  MR  relaxation  times.  The  in¬ 
herent  parameters:  proton-Ti,  proton-T2  are  delineated  for  the  four  physio-mechanical 
cartilage  states:  normal-uncompressed,  normal-compressed,  proteoglycan-depleted- 
uncompressed,  and  proteoglycan-depleted-compressed. 

Each  parameter  was  evaluated,  using  a  two-tailed  t-test,  to  determine  if  proteogly¬ 
can  depletion  had  a  significant  effect  on  any  of  the  MR  time  constants.  The  resulting  P 
values  are  given  in  the  last  column  of  table  4.1. 

Since  the  relaxation  time  constants  throughout  the  imaging  slice  of  the  cartilage  speci¬ 
mens  were  not  uniform  prior  to  compression,  we  have  calculated  Ti  and  T2  histograms  of 
the  normal  and  degraded  cartilage  specimens  prior  to  and  after  compression.  Ti  histo¬ 
grams  for  both  specimens  before  and  after  compression  are  shown  in  figures  4.8-4.9.  We 
should  note  that  experimental  factors  could  have  contributed  to  the  relaxation  time  distri¬ 
butions.  For  example,  the  width  of  the  Ti  histogram  could  represent  noise  in  the  Ti  map 
rather  than  represent  an  actual  Ti  distribution,  especially  since  Ti  in  cartilage  has  been 
shown  to  be  relatively  homogeneous  (Xia  et  al.  1998).  When  testing  this  notion  by  evalu¬ 
ating  the  pixel  by  pixel  reproducibility  of  Ti  maps  of  cartilage  specimens  and 
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There  is  no  evidence  that  proton-Ti  changes  with  proteoglycan-depletion,  and  scarce  evi¬ 
dence,  if  any,  that  proton-T2  changes. 

Compressing  cartilage  provided  a  means  for  introducing  a  proteoglycan-depletion 
induced  change  in  proton-T2.  In  the  compressed  state  (using  the  same  pressure),  the  pro- 
ton-T2  of  the  proteoglycan-depleted  cartilage  is  shorter  than  that  of  normal  cartilage.  Fur¬ 
thermore,  the  change  in  proton-T2  associated  with  compression  is  substantially  more  pro¬ 
nounced  for  proteoglycan-depleted  cartilage  than  it  is  for  normal  cartilage. 

The  data  also  affords  an  evaluation  of  the  changes  in  proton  relaxation  time  con¬ 
stants  associated  with  compression.  To  compliment  the  values  displayed  in  table  4.1,  we 
report  the  statistical  significance  of  the  compression  associated  changes  in  relaxation  time 
in  table  4.2. 

4.3.2  Relaxation  Maps 

Spin-lattice  (TO  and  spin-spin  (T2)  magnetic  resonance  relaxation  time  constant  maps 
were  generated  from  a  series  of  slice  selective  spin-echo  images  acquired  both  before  and 
after  compression  of  the  cartilage.  These  maps  are  shown  in  figures  4.4-4.7  for  represen¬ 
tative  normal  and  trypsin  degraded  cartilage  specimens.  In  contrast  to  the  proton  density 
weighted  images  shown  in  chapter  3,  the  relaxation  time  maps  had  a  uniform  intensity 
only  after  compression.  Prior  to  compression,  Ti  was  fairly  heterogeneous  throughout  the 
imaging  slice,  and  high  values  of  T2  were  predominantly  located  in  the  region  of  the  ar¬ 
ticular  surface  in  both  the  normal  and  trypsin  degraded  cases,  agarose  phantoms,  we 
found  relatively  little  consistency.  However,  the  fact  that  the  Ti  distribution  narrows  upon 
compression  could  implicate  the  existence  of  Ti  heterogeneity.  Although  there  was  a  re- 
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duction  in  Ti  for  both  specimens  upon  compression,  the  change  in  Ti  for  the  degraded 
specimen,  upon  compression,  was  substantially  more  dramatic  than  that  corresponding  to 
the  normal  specimen. 

T2  histograms  for  both  specimens  before  and  after  compression  are  shown  in  fig¬ 
ure  4.10-4.11.  Unlike  the  case  for  Ti,  The  reduction  in  T2  for  both  the  normal  and  de¬ 
graded  specimens  were  similar.  However,  T2  measurements  for  the  degraded  specimen 
both  before  and  after  compression  were  slightly  higher  than  those  for  the  normal  speci¬ 
men.  As  is  indicated  in  the  images  of  figures  4.5  and  4.7,  the  histograms  of  figures  4.10- 
4.11  likewise  indicate  that  both  relaxation  time  constants  become  more  uniform  upon 
compression  for  normal  and  degraded  cartilage. 

From  the  Ti  and  T2  histograms,  we  have  calculated  average  values  of  Ti  and  T2 
prior  to  and  post  compression  for  both  the  normal  and  degraded  specimens.  These  values 
are  given  in  Table  4.3.  Notice  that  the  values  in  table  4.3  differ  substantially  from  the  cor¬ 
responding  values  of  table  4.1.  The  methods  for  measuring  relaxation  times  are  different 
in  each  case.  Also,  the  values  in  table  4.1  include  contributions  from  the  fluid  contained 
within  the  Pyrex  filter. 

4.4  Discussion 

The  central  theme  of  this  experiment  was  to  identify  proteoglycan  depletion  in¬ 
duced  trends  in  proton  MR  relaxation  times.  These  trends  could  serve  as  potential  con¬ 
trast  mechanisms  that  indicate  early  cartilage  degeneration  and  thus  early  signs  of  os¬ 
teoarthritis.  Since  the  normal  functionality  of  articular  cartilage  involves  compression,  we 

have  measured  the  MR  relaxation  times  of  articular  cartilage  in  both  uncompressed  and 
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mechanically  compressed  states.  One  motivation  for  the  introduction  of  mechanical  com¬ 
pression  is  that  it  has  been  hypothesized  that  cartilage  softening  may  be  the  earliest 
known  indicator  of  osteoarthritis.  The  establishment  of  cartilage  compression  and  pro¬ 
teoglycan  depletion  permitted  an  evaluation  of  proteoglycan  depletion  induced  changes  in 
proton  and  sodium  MR  characteristics  of  cartilage  in  the  normal  state  and  during  me¬ 
chanical  compression. 

The  results  of  the  present  study  have  implications  for  clinical  research  in  cartilage 
imaging.  Standard  clinical  proton  imaging  that  employs  basic  Ti  and  T2  contrast  mecha¬ 
nisms  is  unlikely  to  distinguish  proteoglycan  depletion  in  articular  cartilage  since  we 
have  demonstrated  that  proton  Ti  and  T2  do  not  vary  significantly  with  proteoglycan  de¬ 
pletion.  However,  we  have  demonstrated  that,  when  cartilage  is  in  a  mechanically  com¬ 
pressed  state,  a  significant  change  in  proton  T2  does  occur  with  proteoglycan  depletion. 
This  result  indicates  that  it  may  be  possible  to  identify  regions  of  articular  cartilage  that 
have  a  low  proteoglycan  content  with  T2  weighted  proton  imaging,  if  the  cartilage  of  in¬ 
terest  is  in  a  compressed  state.  A  limitation  worth  mentioning  is  that  compressed  state 
proton  T2  changes  associated  with  proteoglycan  depletion  are  on  the  order  of  magnitude 
of  their  standard  deviation  among  samples.  Therefore,  relaxation  time  constant,  alone, 
may  provide  information  on  relative  proteoglycan  depletion  (i.e.  a  measure  of  cartilage 
degradation  over  time),  but  will  require  additional  information  to  make  an  absolute  pre¬ 
diction  of  proteoglycan  depletion.  In-vivo,  mechanical  compression  of  articular  cartilage 
can  be  induced  physiologically  by  having  subjects  perform  certain  exercises  prior  to  im¬ 
aging.  Experimental  evidence  supports  the  notion  that  the  rate  of  cartilage  recovery  is 
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slow  enough  to  permit  MRI  measurements  of  cartilage  thickness  dynamics  during  recov¬ 
ery  (Eckstein  et  al.  1998). 

In  addition  to  clinical  implications,  the  results  reported  in  this  study  are  relevant 
to  understanding  the  physiological  mechanisms  of  cartilage  deformation  and  NMR  re¬ 
laxation  mechanisms  in  cartilage.  For  example,  the  proton  relaxation  time  constants  re¬ 
ported  here  are  consistent  with  those  of  previous  work  where  cartilage  compression  was 
performed  via  an  osmotic  technique  prior  to  MR  evaluation  at  atmospheric  pressure 
(Lusse  et  al.  1995).  The  consistency  of  the  present  results  with  those  cited  indicates  that 
changes  in  compression-induced  proton  MR  characteristics  can  be  generated  by  both  me¬ 
chanical  or  osmotic  pressure  gradients.  This  conclusion  has  implications  for  in-vitro  pro¬ 
ton  MR  cartilage  studies  in  that  investigators  should  be  careful  to  match  the  osmolarity  of 
cartilage  tissue  culture  media  with  that  of  synovial  fluid.  A  further  implication  is  that  it  is 
possible  that  baseline  measurements  of  cartilage  proton  relaxation  times  in  vivo  can  be 
influenced  on  the  osmolarity  by  a  subject's  synovial  fluid. 

With  the  assumption  that  proteoglycan-depletion  is  the  primary  cause  of  cartilage 
softening,  these  results  indicate  that  it  may  be  possible  to  infer  cartilage  softening  via  an 
indirect  measurement  of  proteoglycan  content,  by  measuring  cartilage  proton  MR  relaxa¬ 
tion  time  constants  during  mechanical  compression. 

The  experiments  described  in  this  chapter  also  expand  upon  earlier  progress  to 
identify  early  manifestations  of  osteoarthritis  using  magnetic  resonance  imaging  such  as 
observations  of  MR  signal  intensity  changes  across  the  depth  of  the  cartilage  (Paul  et  al. 
1993),  and  measurements  of  MR  relaxation  time  constants  of  normal  and  degraded  card- 
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lage  (Lehner  et  al.  1989).  MR  relaxation  times  have  also  previously  been  observed  to 
change  in  a  depth  dependent  fashion  in  normal  articular  cartilage  specimens  during  com¬ 
pression  (Rubenstein  et  al.  1996).  Our  present  observations  of  the  MR  properties  of  bo¬ 
vine  specimens  of  articular  cartilage  are  consistent  with  these  previous  studies.  However, 
we  were  unable  to  observe  consistent  depth  dependent  results  with  all  specimens.  This 
inconsistency  could  result  from  the  likelyhood  that  juvenile  cartilage  specimens  (which 
we  used)  are  likely  to  be  structurally  distinct  from  adult  cartilage  (Torzilli  et  al.  1988).  In 
addition,  with  our  technique,  it  was  difficult  to  assess  the  MR  characteristics  of  the  carti¬ 
lage  at  the  articular  surface  since  some  imaging  artifacts  did  arise  from  the  interface  of 
the  cartilage  with  the  glass  filter. 
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signal  intensity  (a.u.) 


Figure  4.1.  Typical  proton  spin-lattice  relaxation  time  constant  data  fit  to  exponential 
functions. 
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signal  intensity  (a.u.) 


echo  time  (msec) 


Figure  4.2.  Typical  spin-spin  relaxation  time  constant  data,  prior  to  compression,  fit  to 
an  exponential  function. 
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signal  intensity  (a.u.) 


Figure  4.3.  Typical  spin-spin  relaxation  time  constant  data,  post-compression,  fit  to  an 
exponential  function. 
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Full  Compression 

Normal 

PG  depleted 

P 

proton  Ti  (sec) 

N 

1.54  +  0.06 

1.52  +  0.09 

0.658 

proton  TI  (sec) 

Y 

0.99  +  0.10 

1.07  +  0.09 

0.217 

%change 

-35.6  ±  7.9 

-29.5  ±  5.8 

0.205 

proton  T2  (msec) 

N 

82.4  ±  2.7 

90.6+10.6 

0.155 

proton  T2  (msec) 

Y 

38.8  ±  2.6 

:  29.6  ±  6.2 

0.028 

%change 

-52.9  ±  2.5 

-67.1+6.8 

0.007 

Table  4.1.  Proton  relaxation  parameters.  Relaxation  times  are  given  for  proton  Ti,  and 
proton  T2  for  cartilage  in  the  normal  and  mechanically  compressed  state.  A  comparison  is 
drawn  between  these  parameters  for  normal  cartilage  and  cartilage  that  has  been  enzy¬ 
matically  treated  to  induce  proteoglycan  depletion.  In  each  case,  five  cartilage  specimens 
were  used.  The  relaxation  times  reported  are  the  corresponding  means.  In  addition  to  re¬ 
laxation  times,  the  per-cent  change  between  the  normal  and  compressed  state  is  given. 
For  each  normal  vs.  proteoglycan-depleted  comparison,  the  results  of  a  two-tailed  t-test  is 
given  in  the  rightmost  column  labeled  "P";  the  value  in  this  column  represents  the  prob¬ 
ability  that  the  normal  and  proteoglycan-depleted  values  could  have  arisen  from  sampling 
a  single  population. 


Normal/Degraded 

Relaxation  Time  Constant 

P  value 

Normal 

T1 

1.60E-05 

Normal  ' 

12 

5.00E-09 

Degraded 

T1 

4.00E-05 

Degraded 

!  .  Jt2  .  : 

3.20E-05 

Table  4.2.  The  statistical  significance  of  the  change  in  relaxation  time  constants  associ¬ 
ated  with  compression. 
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— 

before 

Figure  4.4  Two-dimensional  slice  selective  Ti  maps  of  cartilage  specimen  before  and 
after  compression. 
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Figure  4.5  Two  dimensional  slice  selective  T2  maps  of  cartilage  specimen  before  and  af¬ 
ter  compression. 
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Figure  4.6  Two-dimensional  slice  selective  Ti  maps  of  proteoglycan  depleted  cartilage 
specimen  before  and  after  compression. 
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Figure  4.7  Two-dimensional  slice  selective  T2  maps  of  proteoglycan  depleted  cartilage 
specimens  before  and  after  compression. 
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Figure  4.8  Ti  histogram  of  normal  cartilage  specimen. 
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pixel  frequency 


deplete  before 
deplete  after 


proteoglycan  depleted  specimen 
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Figure  4.9  Ti  histogram  of  proteoglycan  depleted  cartilage  specimen. 
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pixel  frequency 
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pixel  frequency 


Figure  4.11  T2  histogram  of  proteoglycan  depleted  cartilage  specimen. 
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Tl(sec) 

T2(msec) 

Normal  Before 

0.99 

66.4 

0.62 

t  if ,,  38.4 

li? .  '*j '  ill 

Degraded  Before 

0.91 

71.9 

Table  4.3  Proton  MR  mean  relaxation  time  constants  of  normal  and  proteoglycan  de¬ 
pleted  cartilage  specimens.  Degradation  was  conducted  via  a  proteolytic  regime  dis¬ 
cussed  in  the  experimental  protocol  section. 
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Chapter  5 


The  Cartilage  Uni-Axial  Confined  Compression  Experiment 
Within  a  Magnetic  Resonance  Imaging  Environment 

5.1  Introduction 

The  cartilage  compressions  demonstrated  in  chapter  3  indicated  static  conforma¬ 
tions  of  a  cartilage  specimen  before  and  after  a  continuous  compressive  event.  The  dy¬ 
namics  of  cartilage  conformation  during  the  deformation  process  would  lend  insight  into 
cartilage  rheological  mechanisms.  This  chapter  describes  a  technique  that  employs  MRI 
to  monitor  cartilage  compression  with  a  temporal  resolution  on  the  order  of  seconds. 

5.2  Methods 

One  method  of  measuring  the  mechanical  properties  of  a  viscoelastic  material  is 
with  the  uniaxial  confined  deformation  creep  experiment.  In  this  method,  a  constant  pres¬ 
sure  is  applied  to  the  viscoelastic  material  as  the  deformation  of  the  material  is  measured 
as  a  function  of  time  (Fung  1965).  In  the  context  of  the  experiment  with  cartilage,  a  first 
step  is  to  keep  the  experimental  parameters  one  dimensional  such  that  the  analytical  as¬ 
pects  are  manageable  (Mow  et  al.  1990).  This  is  primarily  possible  with  cartilage  as  a 
consequence  of  the  tissue’s  natural  deformation  mechanism.  Cartilage  can  be  thought  of 
as  a  substance  that,  in  bulk,  is  incompressible  unless  fluid  extrusion  is  permitted.  We  thus 
confined  the  cartilage  discs  laterally,  while  allowing  fluid  to  extrude  from  the  articular 
surface.  In  this  way,  both  the  pressure  and  the  response  is  contained  in  one  dimension,  the 
direction  perpendicular  to  the  articular  surface. 
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The  uniaxial  confined  deformation  creep  experiment  was  chosen  as  a  method  to 
measure  the  viscoelastic  properties  of  the  cartilage  specimens.  This  kind  of  experiment 
was  first  employed  on  cartilage  specimens  by  Mow  and  coworkers  (Mow  et  al.  1984).  In 
this  experiment,  a  pressure  is  applied  to  the  cartilage  in  the  form  of  a  step  function.  Prior 
to  pressure  application,  a  series  of  two-dimensional  90°-180°  spin-echo  images  were  ac¬ 
quired.  Upon  compression,  one  dimensional  90°-180°  spin-echo  image  projections  were 
acquired  every  5  seconds  for  a  period  of  80  minutes.  The  compression  was  driven  by  a 
pressure  of  0.690  MPa  (lOOpsi)  which  was  established  within  the  first  20  seconds  of  the 
one-dimensional-projection-image  acquisition  period,  and  was  maintained  for  the  dura¬ 
tion  of  the  creep  experiment. 

5.3  Results 

To  measure  the  compressive  response  of  the  cartilage  disc  to  the  applied  pressure,  we  de¬ 
rived  the  thickness  of  the  cartilage  specimen  from  the  one-dimensional  90°-180°  spin- 
echo  image  projections  that  were  acquired  every  5  seconds  over  a  period  of  80  minutes. 
Six  of  these  projections  are  shown  in  figure  5.1  as  examples.  The  presence  of  signal  in¬ 
tensity  in  these  projections  indicates  the  presence  of  protons,  and  for  the  most  part  indi¬ 
cates  the  presence  of  the  cartilage.  Toward  the  left  hand  side  of  the  peak  is  a  low  intensity 
plateau  that  corresponds  to  the  water  within  the  pores  of  the  permeable  glass  filter.  The 
representation  of  this  projection  does  not  have  sharp  edges.  Two  experimental  factors 
contribute  to  this  lack  of  sharpness.  First,  the  cartilage  specimen  was  not  exactly  perpen¬ 
dicular  to  the  phase  encoding  direction.  The  tilt  was  slightly  less  than  3  degrees.  Second, 
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a  susceptibility  artifact  existed  at  the  interface  between  the  cartilage  specimen  and  the 
ceramic  piston  as  well  as  between  that  cartilage  specimen  and  the  glass  filter. 

The  thickness  of  the  cartilage  was  thus  estimated  as  the  full  width  of  the  projec¬ 
tion  at  half  its  maximal  value.  The  thickness  calculation  entailed  a  linear  interpolation  of 
signal  intensity  for  regions  between  those  represented  by  pixels.  The  thickness  calcula¬ 
tion  was  performed  independently  for  each  one-dimensional  projection.  The  one¬ 
dimensional  projection  images  were  not  collected  beyond  the  80-minute  data  acquisition 
period;  it  is  likely,  based  on  the  trend  of  the  creep  curves,  that  the  cartilage  specimens 
would  have  compressed  further  given  more  time. 

We  arrived  at  using  the  full  width  at  half  maximum  for  the  cartilage  thickness  by 
starting  with  the  assumption  that  our  signal  acquisition  is  independent  of  position  within 
the  field  of  view.  We  therefore  convoluted  a  hypothetical  rectangular  object  in  frequency 
space  with  various  passive  filter  functions  putatively  representing  the  effect  of  our  meas¬ 
urement  system  (including  rotations  of  the  rectangular  object  that  correspond  to  the  carti¬ 
lage  tilt)  and  found  that  the  signal  intensity  at  the  positions  corresponded  to  the  edges  of 
the  initial  box  function  always  contained  the  value  of  half  the  maximum  height  of  the  re¬ 
sulting  curve. 

Given  that  we  have  deduced  the  cartilage  thickness  at  5  second  intervals  during 
the  deformation  creep  experiment,  we  calculated  strain  at  each  of  these  time  points,  de¬ 
fined  as  the  fractional  change  in  thickness  of  the  cartilage  from  its  initial  non-compressed 
state  as  described  by  equation  5.1 
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ef  =■ 


L,-L„ 


[5.1] 


where  ei  is  strain  at  time  i,  Li  is  cartilage  thickness  at  time  i,  and  Lo  is  the  initial  cartilage 
thickness. 

In  doing  so,  we  obtained  the  deformation  creep  curves  shown  in  figure  5.2,  which 
are  similar  to  those  found  in  previous  purely  mechanical  studies  (Mansour  et  al.l 976; 
Mow  et  al.  1980).  Most  notably,  Mow  et.  al.  observed  a  15%  strain  for  an  applied  pres¬ 
sure  of  0.1  MPa.  The  shape  of  the  creep  curve  generated  in  Mow’s  experiment  is  similar 
to  the  shape  we  obtain,  and  has  a  similar  time  constant.  We  have  fit  these  curves  to  a  gen¬ 
eral  exponential  decay  function  given  by  equation  5.2, 


e  =  m 


(i-J-f)l 

V  m2  )) 


( 


+  m, 


t  J 

1-exp  - — 

V  m5  )) 


+  m. 


[5.2] 


where  e  is  strain,  t  is  time,  and  (ml5  m2,  m3, 1x14,  m5)  are  fit  parameters.  Two  time  con¬ 
stants  were  required  for  an  adequate  fit.  The  least  squares  fit  results  are  given  in  table  5.1. 


5.4  Discussion 

There  are  a  multitude  of  theories  describing  the  mechanical  properties  of  articular 
cartilage  that  include  the  electrostatic  interactions  among  negatively  charged  proteogly¬ 
cans  (Buschmann  et  al.  1995),  the  viscous  interaction  between  the  cartilage  fluid  and  its 
solid  matrix  (Mow  et  al.  1984),  and  the  microstructural  organization  of  cartilage 
(Schwartz  et  al.  1994).  It  has  been  proposed  that  purely  viscoelastic  models  of  articular 
cartilage  may  not  be  adequate  for  a  comprehensive  treatment  of  the  mechanical  properties 
of  cartilage  (McCutchen  et  al.  1982).  However,  the  distinction  between  poroelastic  and 
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viscoelastic  deformations  could  be  a  matter  of  semantics.  Viscoelastic  models  are  found 
to  be  consistent  with  numerous  experimental  observations  (Kovach  et  al.  1996). 

We  chose  to  apply  a  simple  rheological  model  to  our  data  and  expand  it  as  neces¬ 
sary.  Since  the  data  from  our  deformation  creep  experiment  at  first  glance  approximates  a 
function  that  exponentially  rises  to  a  fixed  value,  we  chose  the  Voigt  (Fung  et  al.  1965) 
viscoelastic  model  shown  in  figure  5.3  which  is  described  by  equation  5.3. 


where  e  is  strain,  P  is  pressure,  E  is  the  elastic  modulus,  v  is  the  viscosity,  and  c  is  an  ar¬ 
bitrary  constant.  However,  we  found  that  two  time  constants  were  necessary  for  an  ade¬ 
quate  fit  to  our  observations.  Thus,  the  mathematical  description  of  two  Voigt  solids  in 
series  adequately  fit  the  data  as  shown  in  figure  5.2. 

Our  experiment  was  designed  to  mimic  physiological  pressures  which  are  well 
beyond  the  elastic  limit,  which  does  place  limits  on  the  present  physical  interpretation  of 
a  viscoelastic  model.  However,  we  can  state  that  the  deformation  creep  observations  are 
consistent  with  the  notion  that  the  cartilage  rheology  can  be  described  by  a  linked  defor¬ 
mation  of  two  components. 

In  general  the  observed  changes  in  permeability  and  elastic  moduli  resulting  from 
degradation  are  consistent  with  mechanical  studies  involving  alternative  arthritic  models, 
such  as  a  canine  surgically  induced  arthritic  model  (Setton  et  al.  1994),  as  well  as  other 
enzymatic  degradation  models  (Bader  et  al.  1994).  Furthermore,  our  observations  that 
permeability  increases  with  trypsin  digestion  is  consistent  with  the  notion  that  a  reduction 
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in  fixed  charge  density  within  the  cartilage  matrix  induces  an  increase  in  permeability 
(Guetal.  1993). 
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Figure  5.1  One  dimensional  90°-180°  spin-echo  projection  images  of  the  cartilage  speci¬ 
men.  Six  projections  are  shown  at  various  times  after  onset  of  pressure  application  as  in¬ 
dicated  in  the  figure.  The  plateau  on  the  left  corresponds  to  water  held  within  the  glass 
filter  shown  in  figure  1.  The  left  hand  slope  of  the  projection  corresponds  to  the  articular 
surface  of  the  cartilage  specimen.  The  right  hand  side  of  the  projection  corresponds  to  the 
cartilage-bone  boundary.  The  compressive  force  is  from  right  to  left  against  the  stationary 
glass  filter. 
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Figure  5.2  Deformation  creep  curve  of  normal  and  degraded  cartilage  specimens.  Strain, 
defined  as  the  fractional  change  in  cartilage  thickness  resulting  from  a  constant  pressure 
application  of  0.690  MPa,  is  shown  as  a  function  of  time,  the  solid  lines  represent  biex¬ 
ponential  least  squares  fit  to  the  cartilage  thickness  calculations,  represented  as  circles. 
The  curve  for  the  normal  cartilage  specimen  is  labeled  (n),  and  the  curve  corresponding 
to  the  trypsin  degraded  cartilage  specimen  is  labeled  (d).  The  standard  deviations  of  strain 
measurements  are  shown  as  dotted  lines. 
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1  parameter 

normal 

degraded 

ml 

0.183 

0.275 

29.2 

m3 

-Q31 

0.080 

Hisim 

m5 

3.71 

1.57 

Table  5.1  The  deformation  creep  data  of  figure  5.3  were  fit  to  the  biexponential  decay 
function  given  in  equation  5.2. 
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Figure  5.3  Voigt  viscoelastic  model.  This  basic  viscoelastic  model  consists  of  an  elastic 
spring  arranged  in  parallel  to  a  viscous  dashpot.  As  indicated,  E  is  the  elastic  modulus, 
and  v  is  the  viscosity. 


71 


Chapter  6 


Articular  Cartilage  Fluid  Permeability  During  Compression 

6.1  Introduction 

Second  to  bulk  modulus,  the  most  fundamental  parameter  that  describes  the  de¬ 
formation  characteristics  of  cartilage  is  permeability.  The  ease  with  which  fluid  is  able  to 
travel  through  the  cartilage  matrix  dictates  how  quickly  or  slowly  cartilage  deforms  in 
response  to  an  applied  load.  A  first  approximation  of  the  way  that  cartilage  permeability 
affects  the  cartilage  deformation  rate  is  described  in  the  previous  chapter  by  equation  6.3. 
Cartilage  permeability  can  be  derived  from  the  cartilage  deformation  creep  curves  of 
chapter  6.  This  chapter  describes  the  method  for  doing  so,  and  also  provides  an  example 
calculation  that  corresponds  to  the  creep  curves  of  figure  6.2. 

6.2  Permeability  Calculation 

The  data  from  the  deformation  creep  experiment  permits  us  to  calculate  the  per¬ 
meability  of  the  cartilage  assuming  that  the  deformations  are  within  the  elastic  limit.  The 
most  direct  approach  to  calculating  permeability  within  this  experimental  context  is  to 
use  Darcy's  Law  which  states  that  the  velocity  of  the  fluid  in  a  permeable  solid  matrix  is 
linearly  proportional  to  the  pressure  gradient  (Mansour  et  al.  1976).  The  permeability  of 
the  matrix,  k,  is  defined  as  the  constant  of  proportionality.  This  relationship  is  expressed 
symbolically  in  equation  6.1. 

v  =  kVp  [6.1] 
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where  v  is  the  velocity,  and  p  is  the  pressure.  All  three  variables  in  this  relationship  are 
assumed  to  be  functions  of  both  space  and  time. 


Since  our  experimental  setup  is  designed  to  study  the  deformation  of  cartilage  in  one  di¬ 
mension  (i.e.  the  cartilage  is  confined  such  that  it  cannot  expand  laterally),  we  can  treat 
equation  6.1  in  a  one  dimensional  fashion  where  x  is  position.  Since  our  experiment  col¬ 
lects  data  at  discrete  time  points,  it  is  useful  to  explicitly  include  time  in  a  discrete  fash¬ 
ion  as  indicated  by  equation  6.2  where  the  subscript,  i,  indicates  a  particular  point  in  time. 


k  - 1 — ' — 

'  rdp^ 

Idx  Jj 


[6.2] 


It  will  be  useful  to  express  equation  6.2  completely  in  terms  of  parameters  derived  di¬ 
rectly  from  a  deformation  creep  experiment,  namely:  Lo,  the  original  thickness  of  the  car¬ 
tilage;  P,  the  constant  pressure  applied  to  the  cartilage;  Lf,  the  equilibrium  thickness  of 
the  cartilage  under  pressure  P;  At,  the  time  step  between  each  data  acquisition;  and  e;,  the 
strain  at  time  step  i. 


Since  the  fluid  can  be  considered  incompressible  for  our  purposes,  for  the  sake  of 
simplicity,  we  make  the  approximation  that  the  velocity  of  the  fluid  at  all  points  within 
the  cartilage  is  uniform  and  can  be  expressed  as 


vi  = 


[6.3] 


where 


ALi  =  Li  -  Li-\ 


[6.4] 
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and  Li  is  the  thickness  of  the  cartilage  at  time  i.  Furthermore,  since  our  working  defini¬ 


tion  of  strain  is 


Ls  -L„ 


a  substitution  of  equation  6.5  into  equation  6.3  yields  an  explicit  expression  for  v, 


[6.6] 


where 

=  et  -  e(_!  •  [6.7] 

Let  us  assume  for  simplicity  that  the  elastic  properties  of  the  cartilage  specimen  are  ho¬ 
mogeneous.  With  this  assumption,  and  since  there  are  no  pressure  sources  or  sinks 
within  the  cartilage  volume,  the  pressure  gradient  throughout  the  cartilage  volume  is  con¬ 
stant.  Thus,  the  pressure  can  be  expressed  as  a  linear  function  of  space 

p  =  Ax  +  B  [6.8] 

where  A  and  B  are  defined  by  the  boundary  conditions.  Not  all  of  this  pressure,  however, 
affects  the  velocity  of  the  interstitial  fluid,  since  a  portion,  pe,  is  absorbed  by  the  elastic 
resistance  of  the  cartilage 

Pe  =  Ee  [6.9] 


where  E  is  the  elastic  modulus  of  the  cartilage.  The  change  in  pressure  across  the  thick¬ 
ness  of  the  cartilage  is  thus 

Ap=  P-pe  =  P- Ee  [6.10] 

where  P  is  the  constant  applied  pressure.  The  elastic  modulus,  E,  of  the  cartilage  can  be 

defined  by  the  equilibrium  strain  for  a  given  pressure 
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p 


[6.11] 


(Lf  ~Lo)/Lo 


Substituting  equation  [6.1 1]  into  [6.10]  and  dividing  by  the  thickness  of  the  cartilage  at 
time  i  leads  to  the  uniform  pressure  gradient  given  by 


[6.12] 


Substituting  equation  6.12  and  6.6  into  equation  6.2  we  obtain  an  explicit  expression  for 
permeability  in  experimentally  derived  terms 


k,  = 


[6.13] 


The  biexponential  fits  to  the  deformation  creep  data,  when  substituted  into  equa¬ 
tion  6.13  provide  plots  of  permeability  as  a  function  of  strain  given  in  figure  6.1.  The  fig¬ 
ure  indicates  that,  in  the  uncompressed  state,  the  trypsin  degraded  specimen  is  approxi¬ 
mately  twice  as  permeable  to  the  normal  specimen,  and  that  the  increase  in  permeability 
with  respect  to  strain  is  sharper  than  that  of  the  normal  specimen.  The  permeability  of 
the  specimen  in  the  region  of  high  strain  (greater  than  0.3)  is  unlikely  to  be  accurately 
represented  by  this  method  since  the  fluid  velocity  in  this  range  becomes  small.  Addi¬ 
tionally,  in  the  region  of  high  strain,  our  experimental  conditions  exceed  the  elastic  limit 
for  articular  cartilage  as  reported  by  (Mow  et  al.  1980).  Figure  6.2  shows  the  permeabil¬ 
ity,  calculated  in  the  same  manner,  as  a  function  of  time  for  the  deformation  of  both  the 
normal  and  degraded  specimen.  This  figure  shows  that  during  the  course  of  deformation, 
the  degraded  specimen  is  initially  more  permeable  than  the  normal  specimen,  but  after 

approximately  5  minutes  of  sustained  compression,  the  degraded  specimen  loses  its  per- 
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meability  with  respect  to  the  normal  specimen.  This  phenomenon  is  most  likely  the  result 
of  the  degraded  specimen’s  increased  strain. 
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normal  specimen 


Figure  6.1  Permeability  of  the  normal  and  degraded  cartilage  specimens  as  a  function  of 
strain,  defined  as  the  absolute  value  of  the  fractional  change  in  cartilage  thickness  during 
compression. 
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Figure  6.2  Permeability  of  the  normal  and  degraded  cartilage  specimens  as  a  function  of 
time 
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Chapter  7 


Multi-Nuclear  Study  of  the  Cartilage  Osmotic  Model 


7.1  A  Theory  of  Cartilage  Hydrostatic  Pressure 

The  compressive  resistance  of  articular  cartilage  is  provided,  in  part,  by  an  os¬ 
motic  effect  associated  with  the  increased  sodium  content  as  described  in  chapter  2.  The 
hydrostatic  pressure  of  a  dilute  solution  bounded  by  a  semi-permeable  membrane  is  de¬ 
scribed  by  the  ideal  gas  law  (Fermi  1936) 

P  =  jlRT ,  [7.1] 

where  P  is  pressure,  ji  is  the  difference  in  dissolved  particle  concentration  between  the 
inside  and  outside  of  the  membrane,  R  is  the  ideal  gas  constant,  and  T  is  absolute  tem¬ 
perature. 

The  variable  in  equation  7. 1  can  be  expressed  in  term  of  Donnan  equilibria  us¬ 
ing  the  derivations  provided  in  chapter  2.  Specifically  jj.  can  be  approximated  by, 

p  =  n*-n  +  c*-c,  [7.2] 

where  n*  and  n  are  the  sodium  ion  concentrations  inside  and  outside  of  the  cartilage 
specimen  respectively;  and  c*  and  c  are  the  chloride  ion  concentrations  inside  and  outside 
of  the  cartilage  specimen.  Recall  that  the  normalization  convention  of  chapter  2  (that  the 
bathing  sodium  chloride  concentration  is  one  unit)  indicates 

n  =  c  =  1 ,  [7.3] 

thus 

H  =  n*+c*~ 2.  [7.4] 
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Equations  2.4, 2.5,  and  2.9,  when  combined,  indicate  that 


[7.5] 


where/is  the  cartilage  fixed  charge  density,  and  s  is  the  excluded  volume  fraction.  Figure 
7.1  illustrates  a  numerical  graph  of  equations  7.5  and  7.6. 


The  cartilage  mechanical  response  to  pressure  can  be  derived  as  the  osmotic  con¬ 
sequence  of  the  increase  in  fixed  charge  density  associated  with  a  cartilage  fluid  volume 
reduction.  A  cartilage  specimen  in  the  uncompressed  state  has  an  initial  fixed  charge  den¬ 
sity/in  a  fluid  volume  of  1-5.  As  the  cartilage  specimen  is  longitudinally  compressed  by 
e,  expressed  as  a  fractional  amount  of  the  uncompressed  thickness,  the  cartilage  specimen 
fluid  volume  becomes  1-s-e.  Thus,  modifying  equation  7.6  to  take  into  account  confined 


compression  yields 


Jt  = - x+Jl  +  f  /  — - T 

2(1- j-c)  V  I  2(l - 5 - e) 


A  graphical  illustration  of  equations  7.5  and  7.7  is  given  in  figure  7.2.  Equation  7.7  illus¬ 
trates  the  sodium  and  chloride  ion  balance  of  cartilage  specimens  during  the  confined 
compression  experiment  as  a  function  of  pressure.  This  balance  provides  a  framework  for 
testing  the  cartilage  ionic  model  via  interleaved  sodium  and  proton  MRI. 
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7.2  Methods 


7.2.1  Overall  Protocol 

Cartilage  specimens  were  compressed  in-vitro  with  a  constant  0.69  MPa  pressure 
as  in  chapters  3  through  6.  However,  in  the  present  case,  both  sodium  and  proton  MR 
projection  images  were  acquired  in  an  interleaved  fashion.  MR  relaxation  time  constants 
of  the  two  nuclei  were  measured  prior  to  and  post  compression. 

7.2.2  Hardware  Modifications 

The  same  MR  compatible  pressure  cell,  described  in  detail  in  chapter  3,  was  em¬ 
ployed  for  the  interleaved  sodium/proton  acquisition  experiments.  However,  the  rf  coil 
was  modified.  A  double  tuned  solenoid  rf  coil  was  incorporated  into  the  pressure  cell. 
The  transmit/receive  solenoidal  RF  coil  was  double  tuned  to  the  Larmor  frequencies  of 
(86.1  MHz)  and  23Na  (22.7  MHz)  as  described  by  Schnall  et  al.  1985. 

7.2.3  Proton  and  Sodium  Relaxation  Time  Measurements 

The  standard  inversion  recovery  technique  was  employed  to  obtain  proton  Ti 
measurements.  For  normal  cartilage  specimens,  the  inversion  time  was  varied  with  step 
times  of  200  msec  before  compression  and  100  msec  during  compression.  For  the  pro¬ 
teoglycan  depleted  cartilage  specimens,  the  inversion  time  varied  with  step  times  of  250 
msec  before  compression  and  125  msec  during  compression.  In  each  case  30  free  induc¬ 
tion  decays,  (FIDs)  were  acquired.  Each  FID  contained  1024  data  points  at  a  sampling 
rate  of  0.1  msec. 

To  measure  proton  T2,  the  standard  Hahn  spin-echo  technique  (Hahn  1950)  was 

employed.  For  normal  cartilage,  the  echo  time  (x)  was  varied  by  25  msec  before  com- 
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pression  and  10  msec  during  compression.  For  the  proteoglycan-depleted  cartilage,  x  was 
varied  by  20  msec  before  compression  and  8  msec  during  compression.  Spectra  at  30 
time  points  were  acquired  for  each  T2  measurement.  For  each  time  point,  200  FIDs  were 
averaged  to  obtain  a  spectrum. 

Although  sodium  exhibits  bi-exponential  Ti  in  biological  tissues,  due  to  the  small 
difference  between  the  two  Tis,  it  is  difficult  to  measure  each  independendy.  Therefore, 
measurements  of  sodium  spin  lattice  relaxation  time  constants,  T 1,  were  carried  out  using 
the  same  standard  inversion-recovery  methods  used  for  proton.  However,  since  the  Ti  of 
sodium  is  several  orders  of  magnitude  shorter  than  that  of  proton,  the  inversion  time  was 
varied  with  a  step  time  of  6  msec  before  compression  and  3  msec  during  compression  for 
both  normal  and  proteoglycan-depleted  cartilage  specimens.  Thirty  inversion  time  points 
were  acquired  for  each  Ti  measurement;  200  FIDs  were  acquired  and  averaged  for  each 
inversion  time.  Each  FID  contained  1024  data  points  with  sampling  rate  of  0.1msec. 

The  transverse  relaxation  of  sodium  is  bi-exponential,  and  thus  has  two  time  con¬ 
stants:  T2f  and  T2s.  It  has  been  shown,  using  triple  quantum  filtered  spectra,  that  only  the 
slow  component,  T2s,  is  sensitive  to  compression  (Reddy  et  al.  1996).  Therefore  we  used 
the  same  standard  Hahn  spin-echo  technique  for  the  measurement  of  sodium  transverse 
relaxation  times  as  was  used  as  in  the  case  for  proton  transverse  relaxation  times.  For 
normal  cartilage,  x  was  varied  by  4  msec  before  compression  and  2  msec  during  com¬ 
pression.  For  proteoglycan-depleted  cartilage,  x  was  varied  by  3  msec  before  compres¬ 
sion  and  1  msec  during  compression.  Spectra  at  30  time  points  were  acquired  for  each  T2 
measurement.  For  each  time  point,  200  FIDs  were  averaged  to  obtain  a  spectrum. 
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7.2.4  Projection  Image  Acquisition 

One-dimensional  proton  and  sodium  projection  images,  in  a  direction  perpendicu¬ 
lar  to  the  articular  surface,  were  acquired  with  the  pulsed  gradient  spin  echo  pulse  se¬ 
quence  (PGSE).  The  proton  imaging  parameters  were  TR  =  1  sec,  TE  =  20  msec,  NEX  = 
1.  For  sodium  projections,  the  imaging  parameters  were  TR  =  100  msec,  and  TE  =  3 
msec,  NEX  =  800.  The  acquisitions  of  proton  and  sodium  projection  images  were  inter¬ 
leaved.  The  total  acquisition  time  for  one  set  of  interleaved  proton  and  sodium  projection 
images  was  81  seconds  (one  second  for  the  proton  projection  image  and  80  seconds  for 
the  sodium  projection  image).  100  interleaved  projection  images  were  collected  for  a  to¬ 
tal  duration  of  approximately  135  minutes  during  cartilage  compression. 

7.3  Cartilage  Sodium  Content  Mapping 

7. 3. 1  Effects  of  Compression  on  MR  Relaxation  T imes 

Table  7.1  presents  a  summary  of  the  data  obtained  regarding  the  effects  of  com¬ 
pression  and  PG  depletion  on  cartilage  proton  (from  chapter  4)  and  sodium  MR  relaxa¬ 
tion  times.  The  four  inherent  parameters:  proton  Ti,  proton  T2,  sodium  Ti,  and  sodium  T2 
are  delineated  for  the  four  physio-mechanical  cartilage  states:  normal-uncompressed, 
normal-compressed,  PG-depleted-uncompressed,  and  PG-depleted-compressed. 

Each  of  the  four  parameters  was  evaluated,  using  a  two-tailed  t-test,  to  determine  if  PG 
depletion  had  a  significant  effect  on  any  of  the  MR  time  constants.  The  resulting  P-values 
are  given  in  the  last  column  of  table  7.1.  In  the  uncompressed  state,  the  sodium  MR  time 
constants  show  a  much  stronger  PG-depletion  induced  distinction  than  do  the  proton  MR 
time  constants.  There  is  no  evidence  that  proton-Ti  changes  with  PG-depletion,  and 
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scarce  evidence,  if  any,  that  proton-T2  changes.  On  the  other  hand,  both  sodium-Ti  and 
sodium-T2  increase  with  PG-depletion  to  a  statistically  significant  extent. 

Compressing  cartilage  provides  a  means  for  introducing  a  PG-depletion  induced 
change  in  proton-T2.  In  the  compressed  state  (using  the  same  pressure),  the  proton-T2  of 
the  PG-depleted  cartilage  is  shorter  than  that  of  normal  cartilage.  Furthermore,  the  change 
in  proton-T2  associated  with  compression  is  substantially  more  pronounced  for  PG- 
depleted  cartilage  than  it  is  for  normal  cartilage. 

The  contribution  of  compression  in  distinguishing  the  effects  of  PG-depletion  on 
sodium  MR  relaxation  parameters  is  less  pronounced;  the  change  in  sodium-T2,  associ¬ 
ated  with  compression,  is  similar  for  the  normal  and  PG-depleted  specimens. 

The  data  also  affords  an  evaluation  of  the  changes  in  proton  and  sodium  relaxa¬ 
tion  time  constants  associated  with  compression.  To  compliment  the  values  displayed  in 
table  7.1,  the  statistical  significance  of  the  compression  associated  changes  in  relaxation 
time  are  reported  in  table  7.2. 

7.3.2  Effects  of  Compression  on  Sodium  Concentration 

Several  one-dimensional  sodium  an  proton  projection  images  are  shown  in  figures  7.3 
and  7.4.  The  intensity  of  the  images  can  be  made  approximately  proportional  to  sodium 
ion  and  water  volumetric  concentrations  by  correcting  for  the  effects  of  compression  on 
MR  relaxation  times.  This  approximation  is  conducted  by  interpolating  the  MR  relaxation 
characteristics  during  compression  from  the  beginning  and  end  of  compression  by  con¬ 
sidering  relaxation  times  to  be  proportional  to  cartilage  strain.  This  assumption  has  been 
validated  to  a  reasonable  extent  in  appendix  B  for  proton  MR  relaxation  times,  and  re- 
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mains  speculative  for  the  sodium  ions.  The  image  intensities  are  thus  corrected  for  com¬ 
pression-induced  relaxation  using  equation  7.8. 


[7,8] 

where  /’  is  the  corrected  image  intensity,  I  is  the  original  image  intensity,  Rj  is  the  spin- 
lattice  relaxation  rate,  R2  is  the  spin-spin  relaxation  rate,  TR  is  the  repetition  time,  and  TE 
is  the  echo  time. 

Integrating  the  proton  and  sodium  one-dimensional  image  projection  data  (cor¬ 
rected  via  equation  7.8)  for  each  projection  image  (examples  of  which  are  given  in  fig¬ 
ures  7.3  and  7.4)  provides  a  calculation  of  sodium  ion  concentration  in  the  cartilage  fluid 
component.  The  result  of  this  calculation  is  illustrated  in  figure  7.5. 

Although  figure  7.5  does  not  provide  absolute  sodium  ion  concentrations,  the 
relative  increase  in  sodium  ion  concentration  for  a  given  strain  can  be  evaluated  in  terms 
of  the  theoretical  predictions  described  in  section  7. 1 .  The  information  necessary  for  this 
evaluation  is  the  initial  cartilage  fixed  charge  density.  Spectrophotometric  assays  of  carti¬ 
lage  fixed  charge  density,  as  described  in  the  methods  section  of  chapter  2,  indicate  that 
the  cartilage  specimens  typically  have  fixed  charge  densities  of  approximately  200  mM. 
The  expected  sodium  concentration  can  thus  be  derived  from  the  multiplicative  inverse  of 
equation  7.7,  which  indicates  that,  for  a  fixed  charge  density  of  200  mM,  the  initial  carti¬ 
lage  sodium  concentration  within  the  fluid  compartment  is  291  mM.  This  sodium  ion 
concentration  is  consistent  with  bovine  patellar  cartilage  sodium  concentration  reported 
in  the  literature  (Shapiro  et  al.  2000).  The  multiplicative  inverse  of  equation  7.7  (which 
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provides  the  cartilage  sodium  interior/exterior  sodium  partition)  under  such  conditions  is 
illustrated  in  figure  7.5.  The  equilibrium  predictions  in  figure  7.5  systematically  deviate 
from  the  experimental  results  in  that  the  experimental  cartilage  sodium  concentration 
measurements  are  lower  than  the  equilibrium  condition  during  cartilage  compression. 

7.4  Discussion 

7.4.1  Multinuclear  MR  Relaxation  Characteristics  of  Cartilage  Compression 

As  discussed  in  chapter  4,  standard  clinical  proton  imaging  that  employs  basic  Ti 
and  T2  contrast  mechanisms  is  unlikely  to  distinguish  proteoglycan  depletion  in  articular 
cartilage. 

The  results  presented  in  this  chapter  support  the  notion  that  sodium  imaging  could 
have  a  substantial  impact  on  cartilage  evaluation  in  a  clinical  setting.  We  have  demon¬ 
strated  that  both  sodium  Ti  and  T2  change  significantly  as  a  function  of  proteoglycan  de¬ 
pletion  whether  or  not  the  cartilage  is  in  a  compressed  state.  One  caveat  to  this  result  is 
that  sodium  Ti  and  T2  both  increase  with  proteoglycan  depletion,  and  thus  each  may  off¬ 
set  the  effect  of  the  other  depending  on  the  specifics  of  an  imaging  pulse  sequence.  The 
sodium  relaxation  parameters  that  we  report  can  aid  in  the  development  of  an  optimal  so¬ 
dium  pulse  sequence  for  imaging  cartilage  proteoglycan  loss. 

Using  sodium  MR,  one  can  indirectly  infer  changes  in  PG  content,  without  resort¬ 
ing  to  compression.  Two  limitations  of  using  sodium  MRI  are  that  it  has  a  lower  signal  to 
noise  ratio  and  that  it  requires  stronger  imaging  gradients  to  resolve  thin  objects  such  as 
cartilage.  However,  the  shorter  Tj  of  sodium  can  be  exploited  in  signal  averaging  to  im¬ 
prove  the  signal  to  noise  ratio.  High  quality  sodium  images  of  articular  cartilage  have 


86 


been  acquired  in  vivo  in  a  research  setting  (Reddy  et  al.  1998),  which  lends  promise  that 
the  technology  can  be  moved  toward  clinical  applications.  This  inference  can  then  be  re¬ 
lated  to  cartilage  softening  and  in  turn  to  a  non-compressive  assessment  of  cartilage  func¬ 
tionality. 

7.4.2  Evaluation  of  the  Osmotic  Model’s  Mechanical  Prediction 

The  strain  vs.  time  curve  for  the  cartilage  specimen  used  to  generate  figure  7.5  is 
shown  in  figure  7.6.  This  plot  was  generated  in  the  same  manner  as  was  done  in  chapter 
5.  The  equilibrium  strain,  in  response  to  the  constant  applied  0.69  MPa  pressure,  is  ap¬ 
proximately  0.35.  According  to  the  osmotic  model  developed  here  and  illustrated  in  fig¬ 
ure  7.2,  a  cartilage  strain  of  0.35  corresponds  to  a  pressure  of  0.64  MPa  (assuming  a  200 
mM  initial  fixed  charge  density  and  298  °K).  In  this  case,  the  theoretical  prediction  of 
Donnan  equilibria  is  consistent  with  the  experimental  data  within  8%.  This  result  indi¬ 
cates  that  the  osmotic  contribution  is  a  dominant  contributing  factor  to  cartilage  resis¬ 
tance  to  mechanical  compression. 

7.4.3  Evaluation  of  the  Osmotic  Model’s  Ionic  Prediction 

Figure  7.5  indicates  a  systematic  error  between  predictions  of  the  cartilage  ionic 
equilibrium  model  and  experimental  results.  The  most  likely  reason  for  this  deviation  is 

that  rate  of  ion  diffusion  is  slower  than  the  flow  of  ions  associated  with  compressive  ex¬ 
trusion.  This  explanation  implies  that  without  the  diffusion  of  ions  between  the  extruded 
fluid  and  the  cartilage  interior,  the  ion  concentration  of  the  extruded  fluid  would  be  iden¬ 
tical  to  the  ion  concentration  of  fluid  remaining  within  the  cartilage  volume.  However, 
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substantial  ion  diffusion  does  occur;  the  effects  of  this  diffusion  are  what  give  rise  to  the 
positive  slope  of  the  experimental  data  plotted  in  figure  7.5. 

This  hypothesis  can  be  tested  in  several  ways.  The  experiment  can  be  conducted 
as  a  function  of  temperature.  An  increase  or  decrease  in  temperature  would  have  a  pre¬ 
dictable  effect  on  the  rate  of  ion  diffusion.  The  experimental  equipment  can  be  modified 
to  increase  the  diffusive  communication  between  the  cartilage  and  the  extruded  fluid.  It  is 
plausible  that  the  onset  of  inconsistency  between  the  theoretical  prediction  and  experi¬ 
mental  results  that  occurs  at  a  strain  of  0.15,  results  from  the  onset  of  a  lack  of  interaction 
between  the  cartilage  specimen  and  its  extruded  fluid.  Also,  the  experiment  can  be  ex¬ 
tended  in  time  to  observe  whether,  over  an  extended  time  period,  the  sodium  concentra¬ 
tion  within  the  cartilage  increases  for  a  period  of  time  after  mechanical  equilibrium  is 
reached. 

7.4.4  Cartilage  Dynamics  During  Compression 

The  ionic  model,  as  presented  here,  addresses  mechanical  equilibrium,  but  not 
dynamics.  Cartilage  mechanical  dynamics  are  empirically  presented  in  chapters  5  and  6. 
In  chapter  6,  where  cartilage  compressive  dynamics  are  converted  into  cartilage  perme¬ 
ability,  a  linear  relationship  between  stress  and  strain  was  assumed.  However,  in  the  pre¬ 
sent  chapter,  the  cartilage  stress  strain  relationship  is  slightly  non-linear. 

The  results  of  chapter  5  indicate  that  a  selective  reduction  in  fixed  charge  density 
substantially  reduces  compression  time.  Given  that  the  cartilage  collagen  matrix  re¬ 
mained  mostly  intact,  it  is  likely  that  fixed  charge  density  dominates  permeability  as  a 

determinant  of  cartilage  mechanical  dynamics. 
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Figure  7.1  Cartilage  osmotic  strength  as  a  function  of  cartilage  fixed  charge  density. 
Cartilage  osmotic  strength,  fl,  is  defined  as  the  difference  between  the  cartilage  interior 
and  exterior  total  solute  concentrations.  As  with  the  convention  described  in  chapter  2,  ji 
and  fixed  charge  density  are  normalized  to  the  bathing  sodium  chloride  concentration  for 
algebraic  simplicity.  Using  equation  7.1,  and  a  bathing  sodium  chloride  concentration  of 
137  mM  at  298  °K,  one  unit  of  ji  corresponds  to  a  pressure  of  0.34  MPa. 
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Figure  7.2  Cartilage  osmotic  strength  as  a  function  of  compressive  strain.  The  osmotic 
strength  of  cartilage,  fl,  as  defined  in  the  caption  of  figure  7.1,  is  plotted  as  a  function  of 
compressive  strain.  This  example  calculation  employs  a  typical  cartilage  uncompressed 
fixed  charge  density  of  200  mM.  As  described  in  the  caption  to  figure  7.1,  jU  is  normal¬ 
ized  to  the  bathing  sodium  chloride  concentration,  and  one  ji  unit  corresponds  to  0.34 
MPa  when  the  bathing  sodium  chloride  concentration  is  137  mM  and  the  bathing  tem¬ 
perature  is  298°K. 
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(%change) 

-38.2 

8.2 

-19.1 

7.0 

0.004 

sodium  T2  (msec) 

N  : 

18.0 

3.2 

25/7 

3.8 

0.008 

sodium  T2  (msec) 

Y 

11.2 

0.8 

15.5 

2.1 

0.009 

(%change) 

-36.6 

8.6 

-39.3 

4.1 

0.548 

Table  7.1  Sodium  and  proton  relaxation  parameters.  Relaxation  times  are  given  for  pro¬ 
ton  Ti,  proton  T2,  sodium  Ti,  and  sodium  T2  for  cartilage  in  the  normal  and  mechanically 
compressed  state.  A  comparison  is  drawn  between  these  parameters  for  normal  cartilage 
and  cartilage  that  has  been  enzymatically  treated  to  induce  proteoglycan  depletion.  Stan¬ 
dard  deviations  of  all  relaxation  times  are  given  in  the  adjacent  columns  labeled  "SD"  In 
each  case,  five  cartilage  specimens  were  used.  The  relaxation  times  reported  are  the  cor¬ 
responding  means.  In  addition  to  relaxation  times,  the  per-cent  change  between  the  nor¬ 
mal  and  compressed  state  is  given.  For  each  normal  vs.  PG-depleted  comparison,  the  re¬ 
sults  of  a  two-tailed  t-test  is  given  in  the  rightmost  column  labeled  "P";  the  value  in  this 
column  represents  the  probability  that  the  normal  and  PG-depleted  values  could  have 
arisen  from  sampling  a  single  population. 


91 


Nuclei 

Normal/degraded 

Relaxation 
time  constant 

P  value 

Proton 

Normal 

Ti 

1.60E-05 

Normal 

t2  v: 

5^G0E-G9 

Degraded 

Ti 

4.00E-05 

Degraded 

t2 

3.20E-05 

Sodium 

Normal 

Ti 

3.50E-05 

y3jfftmnai 

■  "% ; 

0.0056 

Degraded 

Ti 

0.018 

:  • 

Degraded 

0.002  | 

Table  7.2.  The  statistical  significance  of  the  change  in  relaxation  time  constants  associ¬ 
ated  with  compression. 
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Figure  7.3.  One-dimensional  90°-180°  spin  echo  proton  projection  images  of  a  cartilage 
specimen.  Five  projections  are  shown  at  various  times  after  onset  of  pressure  application. 
The  compressive  force  is  from  right  to  left  against  the  stationary  glass  filter.  The  left- 
hand  slope  of  the  projection  corresponds  to  the  articular  surface  of  the  cartilage  specimen. 
The  right  hand  slope  of  the  projection  corresponds  to  the  cartilage-bone  boundary. 


93 


Figure  7.4.  One-dimensional  spin-echo  sodium  projection  images  of  a  cartilage  speci¬ 
men.  Five  projections  are  shown  at  indicated  times  after  the  onset  of  pressure  application. 
The  specimen  orientation  is  identical  to  that  of  figure  7.3. 
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Figure  7.5.  Theoretical  prediction  and  experimental  measurements  of  cartilage  sodium 
content  in  the  fluid  compartment  as  a  function  of  mechanical  compression.  The  solid  line 
indicates  the  theoretical  prediction  described  by  equation  7.7.  The  squares  represent  ex¬ 
perimental  data  acquired  during  the  interleaved  sodium/proton  projection  images.  Sodium 
content  is  normalized  to  the  initial  cartilage  sodium  concentration  in  the  fluid  compart¬ 
ment.  Note  that  this  normalization  convention  is  different  than  that  used  for  figures  7.1 
and  7.2. 
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Appendix  A 


Magnetic  Resonance  Imaging  of  Unconfined  Cartilage  Compression 

A.1  Introduction 

The  compression  experiments  of  the  earlier  chapters  describe  cartilage  compres¬ 
sion  as  a  unidirectional  event.  In  these  in-vitro  experiments,  the  cartilage  specimens  were 
confined  laterally  such  that  no  cartilage  expansion  could  occur  during  compression.  This 
experimental  constraint  was  introduced  to  simplify  the  deformation  analysis  by  restricting 
stress  and  strain  to  one  dimension.  However,  physiological  cartilage  deformations  are  not 
likely  to  behave  in  this  manner.  It  is  therefore  also  useful  to  study  cartilage  compression 
during  scenarios  where  cartilage  specimens  are  not  laterally  constrained. 

A.2  Methods 

Bovine  articular  cartilage  specimens  were  obtained  from  bovine  patella  as  de¬ 
scribed  in  section  2.4.2.  These  cartilage  disks  were  cut  with  a  diameter  of  5  mm.  Carti¬ 
lage  compression  was  performed  as  a  step  function  in  time,  with  a  pressure  of  0.69  MPa. 
The  pressure  cell  was  the  same,  and  thus,  the  diameter  of  each  cartilage  specimens  was 
smaller  than  the  inner  diameter  of  the  pressure  cell. 

Images  were  obtained  in  a  4  Tesla  field,  using  a  commercially  available  clinical 
spectrometer  (GE  Signa).  The  cartilage  specimens  were  oriented  in  the  cartilage  com¬ 
pression  device  such  that  the  cartilage  articular  surface  was  parallel  to  the  main  magnetic 
field.  Body  gradients  were  employed  for  imaging;  the  field  of  view  was  4.0  cm.  Slice  se¬ 
lective  two-dimensional  90°-180°  spin  echo  images  were  acquired  along  a  plane  perpen- 
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dicular  to  the  axis  of  the  main  magnetic  field.  Slice  thickness  was  2.0  mm  and  oriented 
along  a  diameter  of  the  cartilage  disk.  Images  were  collected  with  a  192  by  256  acquisi¬ 
tion  matrix.  Slice  selective  images  were  acquired  prior  to  cartilage  specimen  compres¬ 
sion,  and  at  the  following  time  points  after  compression  was  applied:  1  min,  2  min,  3  min, 
4  min,  5  min,  6  min,  7  min,  8  min,  9  min,  10  min,  11  min,  12  min,  13  min,  18  min,  23 
min,  28  min,  33  min,  38  min,  and  43  min.  Frequency  encoding  was  in  the  direction 
perpendicular  to  the  cartilage  specimen  articular  surface.  The  spin  echo  images  employed 
an  echo  time  of  20  msec,  and  a  repetition  time  of  300  msec. 

The  experimental  protocol  was  repeated  for  a  proteoglycan  depleted  articular  car¬ 
tilage  specimen.  This  specimen  was  obtained  in  the  same  manner  as  discussed  above,  ex¬ 
cept  that  it  was  incubated  with  1 .0  mg/mL  bovine  pancreatic  trypsin  for  three  hours  on 
ice  then  for  one  hour  at  37°C.  The  imaging  protocol  was  the  same  as  above,  except  that 
the  acquisition  matrix  was  256  by  256  and  images  were  obtained  at  the  following  time 
points:  0.0  min,  1.4  min,  2.8  min,  4.5  min,  5.9  min,  7.4  min,  8.9  min,  10.3  min,  11.8  min, 
13.2  min,  14.5  min,  16.0  min,  17.5  min,  18.9  min,  20.4  min,  21.8  min,  23.1  min,  24.5 
min,  26.0  min,  27.4  min,  28.9  min,  30.3  min,  31.7  min,  33.0  min,  and  34.5  min. 

A.3  Results 

Figure  A.l  shows  the  two-dimensional  slice  selective  images  of  a  representative 
cartilage  specimen  prior  to  and  during  compression.  Cartilage  compression  over  time  is 
demonstrated.  The  serial  images  indicate  that  the  cartilage  specimen  bulges  out  slightly  in 
the  radial  direction  during  compression. 
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The  simplest  analysis  of  this  compressive  event  is  two  dimensional,  represented 
by  the  direction  perpendicular  to  the  cartilage  articular  surface  (longitudinal)  and  the  di¬ 
rection  radial  to  the  cylindrical  cartilage  disk  (radial).  Longitudinal  strain  as  a  function  of 
time  was  calculated  by  computing  the  linear  interpolated  full  width  at  half  maxima  of  se¬ 
rial  one-dimensional  profile  images  generated  by  summing  the  image  arrays  along  the 
horizontal  direction.  The  result  of  this  analysis,  for  the  present  case,  is  illustrated  in  figure 
A.2.  As  in  chapter  6  (figure  6.2),  the  longitudinal  strain  data  did  not  fit  a  single  time  con¬ 
stant  exponential  decay  well,  but  did  adequately  fit  a  two  time  constant  exponential  de¬ 
cay.  The  longitudinal  strain  in  the  unconfmed  compression  experiment  for  a  normal  carti¬ 
lage  specimen  was  similar  that  of  a  proteoglycan  depleted  cartilage  specimen  in  a  con¬ 
fined  compression.  This  is  apparent  by  a  comparison  of  figures  6.2  and  A.2. 

Radial  strain  as  a  function  of  time  was  calculated  in  the  same  manner  as  longitu¬ 
dinal  strain  with  the  distinction  that  the  image  arrays  were  summed  along  the  vertical  di¬ 
rection.  The  result  of  this  calculation  is  illustrated  in  figure  A.3  where  radial  strain  is 
plotted  as  a  function  of  compression  duration.  Also  shown  in  figure  A.3  is  the  least 
squares  best  fit  of  a  single  time  constant  exponential  decay.  By  comparison  of  figure  A.3 
with  A.2,  the  relative  rates  of  radial  expansion  is  demonstrably  higher  than  that  of  longi¬ 
tudinal  compression.  Cartilage  softening  as  a  result  of  proteoglycan  depletion  is  demon¬ 
strated  in  figure  A.4.  Figure  A.5  indicates  an  equilibrium  strain  of  approximately  0.7, 
which  is  substantially  higher  then  the  equilibrium  strain  reached  for  cartilage  specimens 
when  compressed  with  the  same  pressure  but  in  a  confined  geometry  (see  chapter  6).In 
addition,  by  inspection  of  figures  6.2  and  A.5,  the  compressive  time  constant  is  demon¬ 
strated  to  be  relatively  shorter.  Figure  A.6  shows  radial  strain  as  a  function  of  time  during 
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compression.  The  radial  expansion  shown  in  figure  A.6  is  not  significantly  different  than 
that  shown  for  the  normal  cartilage  specimen  shown  in  figure  A.  3. 

A.4  Discussion 

It  has  been  previously  hypothesized  that,  during  unconfined  compression,  an  ar¬ 
ticular  cartilage  specimen  will  quickly  expand  in  a  radial  direction  while  slowly  compress 
in  the  longitudinal  direction  (Armstrong  et  al.  1984).  Our  experiments  confirm  this  hy¬ 
pothesis  while  quantifying  the  phenomenon.  It  is  worth  noting  that  the  description  of  ra¬ 
dial  expansion  exhibited  by  figure  A.3  is  merely  a  first  approximation  of  the  compressive 
event.  The  description  makes  the  assumption  that  the  radial  expansion  of  an  articular  car¬ 
tilage  disk  during  compression  is  uniform  with  regard  to  depth  from  the  articular  surface. 
However,  inspection  of  the  images  in  figure  A.3  indicates  that  this  is  not  the  case.  The 
cartilage  specimen  tends  to  bulge  out  non-uniformly. 

Neither  confined  or  unconfined  compression  of  articular  cartilage  specimens  ac¬ 
curately  mimic  normal  cartilage  physiology.  We  have  encountered  evidence  that  intact 
articular  cartilage  exists  in  a  state  of  substantial  residual  lateral  strain.  This  evidence  is 
that  when  a  cylindrical  articular  cartilage  specimen  is  removed  from  bone,  the  resulting 
cylindrical  specimen  has  a  diameter  approximately  5%  larger  that  that  of  the  core  drill. 
Articular  cartilage  thus  has  a  natural  propensity  to  rapidly  expand  upon  removal  from 
bone.  Since  there  is  no  water  lost  during  this  expansion,  the  cartilage  volume  must  remain 
constant  and  thus  the  cartilage  specimens  must  contract  in  the  longitudinal  direction  upon 
removal  from  bone.  This  aspect  of  cartilage  physiology  likely  provides  an  additional 
source  of  stiffness  in-vivo. 
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A  comparison  of  the  unconfined  compression  of  normal  and  proteoglycan  de¬ 
pleted  cartilage  specimens  indicates  that  the  longitudinal  compressive  modulus  of  articu¬ 
lar  cartilage  is  compromised  by  proteoglycan  depletion  while  the  associated  radial  expan¬ 
sion  is  relatively  independent  of  proteoglycan  depletion.  This  contrast  indicates  that  the 
initial  radial  expansion  of  articular  cartilage  when  undergoing  a  longitudinal  compressive 
pressure  is  not  mediated  by  an  osmotic  resistance.  It  is  likely  that  the  early  radial  expan¬ 
sion  (and  associated  early  longitudinal  contraction)  results  from  a  quick  rearrangement  of 
the  cartilage  matrix  and  does  not  involve  substantial  fluid  extrusion. 
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Figure  A.l  Serial  slice  selective  spin  echo  images  of  a  cartilage  specimen  during  com¬ 
pression.  Compression  duration  progresses  from  the  top  left  image  across  the  rows  to  the 
bottom  right  image.  Final  compression  duration  is  43  minutes.  The  cartilage  specimen  is 
oriented  in  the  MR  compatible  compression  device  such  that  the  cartilage  articular  sur¬ 
face  is  in  contact  with  the  porous  glass  filter.  Cartilage  compression  is  imparted  by  a  ce¬ 
ramic  piston  acting  from  the  top  of  each  image  pushing  down.  The  porous  glass  filter  is 
present  in  each  image  as  diffuse  region  of  increased  signal  intensity  beneath  each  carti¬ 
lage  image.  Notice  that  the  glass  filter  has  a  larger  diameter  that  does  the  cartilage  speci¬ 
men,  which  indicates  that  the  cartilage  specimen  had  room  to  expand  radially  during 
compression. 
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cartilage  specimen  during  unconfined  com 
constant  fit  to  the  experimental  data. 


time  (min) 


Figure  A.3  Radial  strain  of  the  normal  cartilage  specimen  during  unconfined  compres¬ 
sion.  The  solid  line  represents  a  single  time  constant  fit  to  the  experimental  data.  The  fit 
parameters  are  indicated  in  the  figure. 
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Figure  A.4  Serial  slice  selective  spin  echo  images  of  a  cartilage  specimen  during  com¬ 
pression.  Images  sequence  is  similar  to  that  of  figure  8.1 . 
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Figure  A.5  Longitudinal  cartilage  strain  of  the  proteoglycan  depleted  cartilage  specimen 
during  unconfined  compression.  The  solid  line  represents  a  two  time  constant  fit  to  the 
experimental  data. 
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Figure  A.6  Radial  strain  of  the  proteoglycan  depleted  cartilage  specimen  during  uncon¬ 
fined  compression. 
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Appendix  B 


Cartilage  MRI  Strain  Mapping  During  Compression 

B.l  Introduction 

Purely  mechanical  studies  of  articular  cartilage  compression  typically  involve  de¬ 
scribing  articular  cartilage  as  a  relatively  homogeneous  medium.  During  a  compressive 
event,  a  description  is  usually  given  for  the  strain  of  the  entire  specimen  without  an  ex¬ 
plicit  description  of  which  regions  (or  layers)  exhibit  the  bulk  of  the  compression.  The 
reason  for  this  lack  of  spatial  specificity  is  that  it  is  difficult  to  isolate  which  regions  of  a 
cartilage  specimen  are  experiencing  strain.  One  method  for  overcoming  this  challenge  is 
to  shave  an  articular  cartilage  specimen  in  a  direction  parallel  to  the  cartilage  articular 
surface,  then  perform  mechanical  experiments  on  each  layer  individually.  This  approach 
has  been  realized  experimentally  in  permeability  studies  (Maroudas  et  al.  1976).  There 
are  two  inherent  limitations  to  this  approach:  articular  cartilage  is  a  thin  material  that  can 
only  be  practically  shaved  into  a  limited  finite  number  of  layers,  and  articular  cartilage  is 
rarely  (if  ever)  flat  and  thus  presents  a  challenge  to  physically  slicing  it  into  parallel  lay¬ 
ers. 

The  experiment  presented  in  this  chapter  is  a  novel  way  to  quantitatively  to  ap¬ 
proach  cartilage  compression  with  a  depth  resolved  analysis.  The  technique  employs 
magnetic  resonance  imaging  to  collect  depth  dependant  one-dimensional  images  of  ar¬ 
ticular  cartilage  during  compression  (much  in  the  same  fashion  as  was  described  in  chap¬ 
ter  6).  However,  in  the  present  protocol,  spin  echo  images  were  repeated  collected  with  a 
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range  of  echo  times.  This  data  collection  permitted  a  depth  dependant  calculation  of  the 
cartilage  spin-spin  relaxation  time  (T2).  A  knowledge  of  T2  with  depth  dependant  resolu¬ 
tion  is  then  employed  to  describe  the  relative  depth  dependant  cartilage  solid  content  dis¬ 
tribution  at  each  time  point  during  compression.  This  solid  content  distribution,  in  com¬ 
bination  with  the  notion  that  only  fluid  leaves  cartilage  during  compression,  can  then  pro¬ 
vide  a  mapping  of  solid  content  redistribution. 

B.2  Methods 
B.2.1  Overall  Protocol 

Cylindrical  cartilage  discs  were  acquired  from  the  flat  surface  of  juvenile  bovine 
patellea.  Each  disc  was  placed  into  the  MR  compatible  nonmagnetic  pressure  cell  (pre¬ 
sented  in  chapter  3)  that  compresses  cartilage  specimens  in  one  dimension  against  a  po¬ 
rous  glass  filter.  Each  disc  was  then  compressed  with  a  0.69  MPa  pressure  for  a  period  of 
one  hour.  During  the  compression,  spin-spin  relaxation  time  constant  (T2)  maps  of  the 
cartilage  specimens  were  generated  at  a  frequency  of  one  map  per  minute.  The  T2  maps 
were  one  dimensional  along  an  axis  perpendicular  to  the  cartilage  articular  surface.  This 
protocol  provided  cartilage  T2  measurements  as  a  function  of  depth  from  the  articular  sur¬ 
face  with  a  one  minute  temporal  resolution  during  compression. 

B.2. 2  Cartilage  Specimen  Preparation 

Cartilage  specimens  were  obtained  from  juvenile  bovine  patellae  obtained  within 
six  hours  post  mortem  then  stored  on  ice.  8.0  mm  cartilage  disks  were  removed  from  the 
flat  region  of  the  bone  and  stored  in  phosphate  buffered  saline  (PBS)  at  0°  C. 
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B.2.3  T2  Projection  Mapping 

One-dimensional  T2  projection  mapping  was  conducted  by  acquiring  one¬ 
dimensional  90°-180°  spin-echo  image  projections  along  a  direction  perpendicular  to  the 
articular  surface,  then  fitting  these  projection  images  to  equation  B.l, 

I  =  Ae-kT°+C  [B.l] 

where  I  is  the  image  intensity  of  a  given  pixel,  k  is  I/T2,  Te  is  the  echo  time,  and  C  is  an 
arbitrary  constant.  Imaging  was  conducted  in  a  2  Tesla  field.  The  cartilage  specimen  was 
oriented  such  that  the  articular  surface  was  parallel  to  main  magnetic  field.  A  solenoid  rf 
coil  was  used.  The  rf  pulses  were  not  slice  selective;  thus  the  projection  images  repre¬ 
sented  the  entire  cartilage  specimen  collapsed  onto  a  line  perpendicular  to  the  surface.  No 
phase  encoding  was  employed.  Frequency  encoding  was  perpendicular  to  the  articular 
surface.  The  magnetic  field  gradient  was  supplied  by  a  custom  built  gradient  set  capable 
of  providing  8  gauss/cm. 

One  T2  projection  map  was  acquired  each  minute.  Each  T2  projection  was  calcu¬ 
lated  from  seven  projection  images  having  the  following  echo  times:  20  msec,  40  msec, 
60  msec,  80  msec,  100  msec,  120  msec,  140  msec,  160  msec,  and  180  msec.  A  five  sec¬ 
ond  delay  existed  between  each  projection  image  acquisition. 

B.2.4  Cartilage  Compression 

Cartilage  compression  was  conducted  with  the  non-magnetic  pressure  cell  de¬ 
scribed  in  chapter  3.  The  pressure  cell  employed  a  pneumatic  method  for  pressing  a  ce¬ 
ramic  piston  onto  the  cartilage  which  was  positioned  against  a  porous  glass  filter.  The 
filter  was  connected  to  a  drain  which  permitted  cartilage  extruded  cartilage  fluid  to  exit 
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the  cell  during  compression.  A  continuous  pressure  of  100  psi  (0.69  MPa)  was  employed 
for  a  duration  of  120  min.  During  each  compressive  event,  T2  projection  maps  were  ac¬ 
quired  at  a  frequency  of  one  projection  map  per  minute. 

B.3  Results 

B.3.1  T2  map  calculations 

Onedimensional  apparent  T2  maps  were  calculated  from  the  one  dimensional  90°- 
180°  spin-echo  projection  images.  The  projection  images,  acquired  with  several  echo 
times,  were  fit  to  equation  B.2 

TE 

I  =  AeJl  +  B ,  [B.2] 

where  I  is  image  intensity,  A  is  a  parameter  proportional  to  the  MR  visible  proton  content 
(in  addition  to  being  related  to  the  particular  imaging  conditions),  Te  is  the  echo  time,  and 
B  is  an  offset  associated  with  image  noise.  Figure  B.l  shows  an  example  series  of  projec¬ 
tion  images;  figure  B.2  illustrates  the  associated  T2  map.  Figure  B.3  illustrates  articular 
cartilage  T2  projections  for  the  same  cartilage  specimen  at  two  different  times  during 
compression.  Sixty  similar  T2  projections  per  specimen  were  available  for  analysis  tem¬ 
porally  separated  by  one-minute  intervals. 

B.3. 2  Cartilage  Water  Content  Calculations 

Since  the  elapsed  time  between  projection  image  acquisitions  was  approximately 
five  times  that  of  previously  measured  values  of  cartilage  spin-lattice  relaxation  times 
(TO  (Kaufman  et.  al.,  1999),  it  is  assumed  that  T!  effects  on  signal  intensity  are  negligi¬ 
ble.  The  major  specimen  dependant  contributions  to  signal  intensity  are  thus  proton  den- 
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sity  and  apparent  spin-spin  relaxation  time  (T2)  effects.  It  is  assumed  that  the  major  pro¬ 
ton  signal  is  associated  with  cartilage  water  content  since  the  protons  associated  with  the 
major  cartilage  biopolymers  will  likely  have  short  relaxation  times  which  inhibit  visibil¬ 
ity  under  the  imaging  protocol  used.  Water  content  within  a  multiplicative  factor  was  thus 
obtained  by  rearrangement  of  equation  B.2.  Figure  B.4  illustrates  the  results  of  an  exam¬ 
ple  water  content  calculation  associated  with  the  data  from  figures  B.l  and  B.2.  Relative 
water  content  before  and  after  one  hour  of  a  continuous  mechanical  pressure  of  0.69  MPa 
is  illustrated  in  figure  B.5. 

B.3.3  Cartilage  Longitudinal  Strain  Calculations 

Cartilage  thickness  was  determined  as  the  full  width  at  half  maximum  of  the  one¬ 
dimensional  cartilage  projection  images  as  validated  previously  (Kaufman  et.  al.,  1999). 
Linear  interpolation  between  adjacent  pixels  was  used  routinely  since  the  position  of  in¬ 
tensity  at  half  maximum  often  occurred  between  pixel  boundaries  in  the  one-dimensional 
images.  Cartilage  longitudinal  strain  was  calculated  as  the  change  in  thickness  of  the  car¬ 
tilage  specimen  divided  by  the  cartilage  specimen’s  original  thickness.  A  typical  cartilage 
specimen  strain  vs.  time  curve  for  a  constant  compressive  force  of  0.69  MPa  is  shown  in 
figure  B.6.  The  cartilage  specimen  used  to  generate  figure  B.6  was  the  same  as  in  figures 
B.1-B.5. 

The  boundaries  generated  while  determining  cartilage  longitudinal  strain  were 
employed  when  calculating  the  total  water  volume  of  a  cartilage  specimen  during  com¬ 
pression.  The  total  cartilage  specimen  water  volume  was  calculated  as  the  integral  of  car¬ 
tilage  water  volume  over  the  area  bound  by  the  cartilage  specimen  thickness  boundaries. 
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Linear  interpolation  was  used  when  calculating  the  water  volume  contributions  of  the 
signal  portions  associated  with  the  cartilage  specimen  edges,  which,  in  general,  did  not 
correspond  exactly  to  pixel  boundaries.  Figure  B.7  illustrated  relative  total  cartilage  water 
content  as  a  function  of  time  during  compression.  Figure  B.8  illustrates  relative  total  car¬ 
tilage  water  content  as  a  function  of  longitudinal  strain. 

The  boundaries  used  to  calculate  longitudinal  strain,  and  relative  total  water  con¬ 
tent  were  also  used  to  calculate  aggregate  cartilage  T2;  the  aggregate  value  of  T2  is  de¬ 
fined  here  as  simply  the  average  of  the  T2  values  at  each  pixel  within  the  cartilage 
boundaries.  Figure  B.9  illustrates  aggregate  T2  as  a  function  of  time.  Figure  B.10  illus¬ 
trates  the  relationship  between  aggregate  T2  and  cartilage  longitudinal  strain.  Figure  B.ll 
illustrates  the  relationship  between  aggregate  T2  and  cartilage  volume.  Figure  B.12  illus¬ 
trates  the  relationship  between  aggregate  T2  and  relative  cartilage  water  content. 

B.3.4  Lateral  expansion  of  Articular  Cartilage 

During  the  compressive  event,  the  cartilage  specimen  experiences  a  reduction  in  volume 
by  the  extrusion  of  a  portion  of  its  fluid  content  as  indicated  by  figure  B.7.  The  cartilage 
specimen  also  exhibits  a  radial  expansion  perpendicular  to  the  direction  of  compression. 
Figure  B.12  exhibits  a  framework  of  terms  useful  for  describing  the  compressive  event. 
In  this  figure,  r  is  the  original  radius  of  the  cartilage  specimen,  Ar  is  the  radial  increase 
associated  with  lateral  expansion,  and  e  is  the  strain  in  the  longitudinal  direction  as  de¬ 
fined  as  the  fractional  change  in  cartilage  thickness.  For  algebraic  simplicity,  all  dimen¬ 
sions  are  normalized  to  the  original  thickness  of  the  cartilage  specimen.  The  original  car¬ 
tilage  volume,  Vo,  is  thus 
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VQ=7tr\ 


[B.3] 


and  the  volume  occupied  by  the  cartilage  solid  component,  Vs,  is  thus 

Vs  =  sJtr2  .  [B.4] 

The  total  cartilage  volume,  VT,  as  a  function  of  longitudinal  (e)  and  lateral  (  Ar)  strain  is 
defined  geometrically  as 

VT(e,Ar)=(l-e)rt(r+Ar)2,  [B.5] 

and  thus  water  content,  W,  is  described  by 

W  =  (l-eyt(r  +  Arf -sjrr2  .  [B.6] 


Since  our  experimental  protocol  indirectly  measured  water  content  and  longitudi¬ 
nal  strain  as  a  function  of  time,  it  is  useful  to  rearrange  equation  B.5  and  express  lateral 
strain  as  the  dependant  variable 


A  r  =  . 


\W  +  snr 2 

I  (1-e)* 


r . 


[B.7] 


The  calculations  employed  in  figures  B.6  and  B.7  can  be  substituted  in  the  right 
side  of  equation  B.7  to  generate  a  prediction  of  the  radial  expansion  function  of  a  carti¬ 
lage  specimen  during  a  constant  compressive  event.  With  the  measured  20%  volume 
fraction  of  the  cartilage  solid  component,  cartilage  radial  expansion  is  calculated  as  a 
function  of  compression  time  and  is  illustrated  in  figure  B.13.  Radial  expansion  was 
measured  directly  during  separate  cartilage  compression  experiments  at  4  Tesla,  where 
two-dimensional  images  were  acquired  at  a  one-minute  resolution.  A  comparison  of  the 
cartilage  radial  expansion  prediction  from  the  one-dimensional  projection  images  to  the 

radial  expansion  measurements  via  two-dimensional  imaging  is  shown  in  figure  B.14. 
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The  theoretical  prediction  is  comparable  to  the  experimental  results,  however,  the  predic¬ 
tion  has  a  temporally  broad  depression  centered  at  35  minutes  after  pressure  onset,  which 
is  not  seen  experimentally. 

B.3.5  Cartilage  Strain  Mapping  During  Compression 

The  results  for  longitudinal  compression  and  radial  expansion  as  a  function  of 
time  provide  a  temporal  description  of  cartilage  compression  in  terms  of  overall  specimen 
morphology.  This  description  combined  with  the  water  content  estimations  described  in 
figures  B.4,  B.5,  and  B.7  can  provide  estimates  of  the  cartilage  compressive  behavior  as  a 
function  of  distance  from  the  cartilage  articular  surface.  The  water  content  normalization 
of  figures  B.4,  B.5,  B.7,  and  B.8  is  based  on  equation  B.5,  where  the  initial  cartilage  wa¬ 
ter  content  must  be  consistent  with  the  conditions  of  zero  initial  strain 

W(g  =  0,  Ar  =  0)  =  (l  -  s)jtr2 .  [B.8] 

The  initial  cartilage  measurements  for  the  example  shown  are  that  the  initial  cartilage 
solid  volume  content  was  20%  (s=0.2)  and  that  the  initial  cartilage  diameter  was  4.67 
times  larger  that  its  thickness  (r=2.335).  Given  that  the  cartilage  geometry  predicts  total 
cartilage  volume  as  a  function  of  time,  the  water  content  calculations  can  be  converted 
into  a  description  of  cartilage  solid  content  flow.  The  solid  content  flow  can  then  be  used 
to  describe  longitudinal  strain  mapping  of  the  cartilage  specimen  during  compression. 
Each  pixel  in  the  one-dimensional  projection  description  of  water  content  corresponds  to 
a  virtual  cylindrical  slice  of  the  cartilage  specimen.  If  the  cartilage  thickness  is  initially 
divided  into  N  pixels,  the  volume  associated  with  a  given  pixel  is 
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VT  = 


_  n{r  +  Ar)2 
~N 


[B.9] 


where  vT  is  the  total  volume  represented  by  a  given  pixel.  Substitution  of  equation  B.7 
into  equation  B.9  yields  an  explicit  expression  for  pixel  volume  as  a  function  of  longitu¬ 
dinal  strain  and  total  cartilage  water  content 


vT  —  ■ 


W  +  sftr 


[B.10] 


N(l  -  e)  ' 

The  corresponding  water  content  for  each  slice  can  be  calculated  from  the  water 
content  maps,  given  that  the  composite  water  content  map  containing  all  time  points  is 
normalized  to  be  consistent  with  equation  B.6  and  B.8.  The  water  content  maps,  thus 
normalized,  when  subtracted  from  the  total  volume  maps  predicted  by  equation  B.10  on 
an  individual  pixel  basis  provide  maps  of  solid  content  for  each  time  point.  The  solid  con¬ 
tent  map  for  the  same  example  specimen  is  shown  in  figure  B.14. 


An  understanding  of  the  movement  of  solid  content  provides  a  way  to  determine 
the  compressive  behavior  of  individual  cartilage  regions.  Computational  demarcations 
were  determined  in  a  direction  perpendicular  to  the  cartilage  articular  surface.  The  solid 
content  corresponding  to  the  initial  cartilage  one-dimensional  projection  image  was  di¬ 
vided  into  ten  contiguous  evenly  distributed  regions.  Total  solid  content  corresponding  to 
each  individual  region  was  calculating  by  integrating  the  initial  (t  =  0)  solid  content  map 
individually  over  the  demarcated  spatial  regions.  Thus,  the  discrete  total  solid  content  per 
region  function  was  determined.  This  process  can  be  expressed  symbolically  by  figure 
B.15  and  equationsB.il -B.14. 
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a  =  lx— 1) — 

v  '10 

[B.12] 

i  N 

b  =  x - 1 

[B.13] 

10 

xe  {1,2,..., 10} 

[B.14] 

where  O(x)  is  total  solid  content  for  each  demarcated  unit,  x  is  the  demarcated  unit  index, 
<j)  is  the  solid  content  for  each  image  pixel,  a  and  b  are  the  pixel  limits  for  each  demar¬ 
cated  unit,  and  N  is  the  total  number  of  pixels  across  the  cartilage  specimen. 

The  discrete  function  <E>(x)  provides  a  means  for  strain  mapping  given  the  follow¬ 
ing  three  assumptions:  cartilage  volume  is  primarily  decreased  during  compression  by  a 
loss  of  fluid  content,  cartilage  retains  a  constant  total  solid  content  during  compression, 
the  longitudinal  ordering  of  cartilage  solid  content  remains  constant  during  compression. 
In  other  words,  it  is  assumed  that  temporal  transformations  of  cartilage  solid  content 
functions  are  homeomorphic.  This  assumption  indicates  the  algorithm  for  mapping  carti¬ 
lage  longitudinal  strain  shown  in  figure  B.16,  where  P(x)  is  the  pixel  position  function  of 
the  cartilage  computational  demarcations.  This  algorithm,  when  applied  to  all  time  points 
during  compression,  provides  P(x,t),  a  description  of  cartilage  position  movements  as  a 
function  of  time. 

P(x,t)  can,  in  turn,  be  converted  into  a  spatial  mapping  of  cartilage  longitudinal 
strain  during  compression 
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P(x,t)-  P(x-\,t) 


[B.15] 


P(x,0)-P(x-l, 0) 

where  ex(t)  indicates  the  longitudinal  strain  of  the  xth  computational  demarcation  (region 
x)  at  time  t.  A  representative  ex(t)  is  shown  in  figure  B.17.  Each  function  approximates  an 
exponential  decay.  Longitudinal  compressive  modulus  and  compressive  time  constant 
can  be  obtained  by  fitting  each  ex(t)  to  equation  B.  16 


.(0=  E, 


l-exd 


f-t" 


[B.  16] 


where  Ex  is  the  compressive  modulus  of  region  x,  and  Tx  is  the  compressive  time  constant 
of  region  x.  The  results  of  fitting  the  experimental  results  to  each  ex(t)  are  illustrated  in 
figures  B.18  and  B.19. 


B.4  Discussion 

In  this  experiment,  the  compressive  properties  of  an  articular  cartilage  specimen 
have  been  deduced  as  a  function  of  distance  from  the  articular  surface.  The  results  indi¬ 
cate  the  a  relatively  soft  and  compliant  region  of  cartilage  exists  between  the  highly  calci¬ 
fied  region  of  cartilage  adjacent  to  bone  and  the  cartilage  central  region.  That  a  soft  re¬ 
gion  of  cartilage  lies  between  the  bone  interface  and  the  main  cartilage  body  makes  sense 
in  retrospect,  but  was  an  unexpected  result. 

Traditionally,  the  softness  distribution  of  cartilage  has  been  described  as  a  gradual 
increase  in  compressive  modulus  in  a  direction  from  the  articular  surface  to  the  cartilage 
bone  interface.  This  description  has  been  supported  by  an  experiment  where  microscopic 
imaging  of  chondrocyte  nuclei  was  conducted  during  compression  (Schinagl  et  al.  1996). 
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The  primary  distinction  between  the  experiment  conducted  by  Schinagl  et  al.  and  the  ex¬ 
periment  presented  in  this  chapter  is  that  of  lateral  expansion.  If  an  experimental  protocol 
does  not  permit  fluid  extrusion  in  a  radial  direction,  then  all  cartilage  fluid  must  exit  the 
cartilage  via  the  articular  surface.  In  other  words,  if  cartilage  fluid  is  retained  laterally,  for 
fluid  to  leave  the  cartilage  region  near  the  bone  interface,  it  must  push  its  way  through  the 
entire  cartilage  depth.  In  this  way,  the  cartilage  specimen  as  a  whole  can  be  no  more  per¬ 
meable  to  fluid  than  the  permeability  of  the  articular  surface. 

The  experimental  protocol  of  this  chapter  permits  cartilage  fluid  to  exit  the  carti¬ 
lage  specimen  at  the  lateral  surface  as  well  as  permits  cartilage  radial  expansion  during 
longitudinal  compression.  The  results  thus  indicate  that,  during  unconfined  compression, 
cartilage  fluid  preferentially  leaves  the  specimen  laterally.  It  is  difficult  to  interpret  the 
implications  of  this  behavior  since,  as  discussed  in  chapter  8,  cartilage  is  likely  to  exist 
in-vivo  with  a  substantial  negative  radial  residual  strain.  Furthermore,  in  a  in-vivo  situa¬ 
tion,  cartilage  fluid  ultimately  must  exit  some  portion  of  the  articular  surface  during  a 
compressive  event  that  incurs  a  total  cartilage  volume  reduction. 

The  strength  of  the  technique  presented  here  is  that  the  present  study  can  be  ex¬ 
tended  to  an  in-vivo  setting.  Novel  methods  for  fast  three-dimensional  T2  mapping  in 
conjunction  with  exercise  induced  in-vivo  cartilage  compression/expansion  experiments 
can  be  combined  with  the  present  analysis  to  obtain  a  true  three-dimension  continuum 
description  of  cartilage  compression  in-vivo. 
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Figure  B.l.  One-dimensional  example  projection  images.  The  direction  of  the  projection 
images  is  such  that  pixel  0  corresponds  to  the  direction  of  the  cartilage  articular  surface 
and  pixel  60  corresponds  to  the  direction  of  the  cartilage  bone  interface.  These  projection 
images  correspond  to  the  cartilage  specimen  several  minutes  after  compression  onset. 
Each  projection  is  of  a  90°-180°  spin-echo  pulse  sequence  with  a  range  of  echo  times  as 
indicated  in  the  figure. 
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Figure  B.3.  T2  before  and  after  compression.  Shown  are  two  one-dimensional  T2  projec¬ 
tions.  The  projections  represent  two  limiting  position  of  the  cartilage  during  compression: 
the  cartilage  T2  projection  prior  to  compression,  and  the  cartilage  T2  projection  after  one 
hour  of  continuously  applied  pressure  of  0.69  MPa. 
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Figure  B.4.  One  dimensional  water  content  map  of  the  cartilage  specimen  after  several 
minutes  of  compression.  The  articular  surface  is  the  left  portion  of  the  map;  the  calcified 
cartilage  region  is  the  right  portion  of  the  map.  the  cartilage  specimen  is  being  com¬ 
pressed  against  a  porous  Pyrex  filter  from  left  to  right. 
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Figure  B.5.  Water  content  before  and  after  compression.  The  water  content  projections 
shown  here  are  similar  to  the  water  content  projection  shown  in  figure  B.4.  These  projec¬ 
tions,  however,  illustrate  the  water  content  states  before  and  after  one  hour  of  a  continu¬ 
ous  0.69  MPa  longitudinal  pressure.  The  time  points  presented  are  the  same  as  those  il¬ 
lustrated  in  figure  B.3. 
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Figure  B.6.  Longitudinal  strain  as  a  function  of  time.  Longitudinal  strain  of  the  cartilage 
specimen  is  presented  as  a  function  of  duration  after  a  continuous  0.69  MPa  longitudinal 
pressure.  The  temporal  resolution  is  one-minute  intervals.  Strain  was  calculated  from  the 
one-dimensional  projection  images  and  is  defined  here  as  the  change  in  cartilage  longitu¬ 
dinal  thickness  divided  by  the  original  cartilage  longitudinal  thickness. 
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Figure  B.7.  Relative  total  water  content  as  a  function  of  time.  Water  content  calculations 
shown  were  calculated  as  integrations  of  the  one  dimensional  water  content  maps  over 
the  cartilage  thickness.  The  integrations  were  normalized  to  maintain  consistency  with 
equation  7. 
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Figure  B.8.  Total  water  content  as  a  function  of  longitudinal  strain.  Water  was  calculated 
as  in  figure  B.7,  and  is  presented  as  a  function  of  cartilage  longitudinal  strain  calculated 
as  described  in  the  legend  of  figure  B.6. 
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Figure  B.10.  Composite  T2  as  a  function  of  longitudinal  strain.  Composite  T2  was  calcu¬ 
lated  as  described  in  the  legend  of  figure  B.9;  longitudinal  strain  was  calculated  as  de¬ 
scribed  in  the  legend  of  figure  B.6. 
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Figure  B.ll.  Composite  T2  as  a  function  of  relative  total  water  content.  Composite  T2 
was  calculated  as  described  in  the  legend  of  figure  B.9;  relative  total  water  content  was 
calculated  as  described  in  the  legend  of  figure  B.7. 
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Figure  B.12.  Side  view  of  articular  cartilage  disk  before  and  during  compression.  The 
solid  line  denotes  the  original  position  of  the  cartilage  disk;  the  dashed  line  denotes  the 
compressed  state.  For  algebraic  simplicity,  the  physical  dimensions  of  the  cartilage  disk 
are  normalized  to  the  original  thickness.  The  original  radius  of  the  cartilage  disk  is  r. 
During  compression,  lateral  (radial)  expansion  occurs  by  Ar.  The  longitudinal  compres¬ 
sion  is  e,  which  as  a  result  of  the  normalization  convention,  is  also  the  longitudinal  strain. 
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Figure  B.13.  Cartilage  specimen  radial  expansion.  The  solid  line  with  square  markers 
represents  the  prediction  of  the  one  dimension  compression  results  using  equation  B.6. 
For  comparison,  the  results  of  a  separate  compression  experiments  are  shown  (dashed 
line)  where  cartilage  lateral  expansion  was  directly  measure  via  two  dimensional  imag- 
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Figure  B.14.  Volumetric  solid  content  of  the  cartilage  specimen.  The  solid  content  was 
calculated  from  the  relationship  (expressed  in  equation  B.9)  between  derivations  of  carti¬ 
lage  volume  and  water  content.  Normalized  was  conducted  to  maintain  consistency  with 
equation  B. 7. 
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Figure  B.15.  Dlustration  of  the  computational  demarcation  scheme.  The  thickness  of  the 
cartilage  specimen  is  divided  into  ten  bins  labeled  {1,2,. ..,10}.  The  longitudinal  pixel 
density  is  N/10  pixels  per  bin. 
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Figure  B.16.  Flow  chart  of  the  algorithm  used  to  calculate  the  locations  of  the  original 
computational  demarcations  at  any  given  point  in  time.  P(x)  is  the  pixel  position  of  de¬ 
marcation  x.  S(x)  is  a  temporary  holding  bin  for  comparing  collected  solid  content  to  the 
original  solid  content  demarcation  array,  d>(x),  introduced  in  equation  B.10. 
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Figure  B.17  Time  dependant  longitudinal  strain  of  cartilage  section.  The  longitudinal 
strain  of  a  representative  cartilage  section  (x  =  1  as  in  figure  B.10)  is  presented  as  a  func¬ 
tion  of  time  under  a  constant  compressive  pressure  of  0.69  MPa.  This  plot  corresponds  to 
d(t)  as  defined  in  equation  B.14.  The  solid  line  is  a  least  squares  best  fit  of  the  derived 
experimental  results  (represented  as  square)  to  equation  B.15.  The  result  of  the  best-fit 
algorithm  is  indicated  in  the  figure.  A  similar  analysis  was  conducted  for  all  sections 

x={  1,...,10},  the  results  of  which  are  displayed  in  figures  B.18  and  B.19. 
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Figure  B.18.  Equilibrium  strain  of  cartilage  sections.  Equilibrium  compressive  strain  is 
plotted  as  a  function  of  cartilage  section  as  demarcated  in  figure  B.15.  The  magnitude  of 
strain  is  proportional  to  compressive  modulus.  The  constant  applied  pressure  is  0.69 
MPa.  The  strain  of  each  section  was  obtained  as  a  best  fit  to  equation  B.15.  Since  strain  is 
assumed  to  be  a  continuous  function  of  cartilage  space,  the  calculated  points  are  con¬ 
nected  by  a  cubic  spline.  Open  squares  indicate  distribution  of  specimens. 
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Figure  B.19  Compression  time  constant  of  cartilage  sections.  Compression  time  constant 
is  plotted  as  a  function  of  cartilage  section  as  demarcated  in  figure  B.15.  The  constant 
applied  pressure  is  0.69  MPa.  The  time  constant  of  each  section  was  obtain  as  a  best  fit  to 
equation  B.15.  Since  compression  time  constant  is  assumed  to  be  a  continuous  function 
of  cartilage  space,  the  calculated  points  are  connected  by  a  cubic  spline.  Open  squares 
indicate  distribution  of  specimens. 
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